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Abstract 

We  have  developed  and  tested  in  the  field  a  deterministic  model  to 
predict  metal  accumulation  in  benthic  organisms.  In  the  formulation  of 
the  model,  we  used  concepts  derived  from  the  free-ion  activity  model 
of  metal-organism  interactions  and  from  surface  complexation  theory, 
to  relate  metal  concentrations  in  the  organisms  to  those  in  the  water 
or  in  the  surficial  oxic  sediments. 

The  model  was  successfully  tested  in  the  field  for  cadmium 
accumulation  in  the  freshwater  bivalve  Anodonta  grandis .    Linear 
regression  analysis  indicated  that  Cd  concentrations  in  the  soft 
tissues  of  the  bivalve,  [Cd(org)],  are  related  to  the   free  cadmium 
concentrations,  [Cd^"*"]  ,  estimated  from  the  total  dissolved  Cd  and 
inorganic  ligands  concentrations:  [Cd(org)]  =  44  [Cd^"*"]  +  10 
(r^=0.81).  In  a  further  step,  partitioning  of  cadmium  between  water 
and  surficial  oxic  sediments  was  interpreted,  using  surface 
complexation  concepts,  in  terms  of  sorption  of  this  metal  on  two 
sediment  components:  organic  matter  and  iron  oxyhydroxides.  This 
partitioning  model  relates  the  free  Cd^"^  concentrations  to  sedimentary 
variables  (e.g. ,  total  sedimentary  Cd,  iron  oxyhydroxides  and  organic 
carbon  concentrations)  and  water  pH;  sedimentary  variables  can  then  be 
used  to  predict  cadmium  concentrations  in  the  tissues  of  A.  grandis. 
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ABIOTIC  FACTORS  INVOLVED  IN 

PREDICTING  TRACE  METAL  LEVELS 

IN  FRESHWATER  BIVALVES 


A.  Tessier  and  P.G.C.  Campbell 


1.  INTRODUCTION. 

Industrialization  has  led  to  increasing  fluxes  of  trace  metals  from  terrestrial  and  atmospheric 
sources  towards  the  aquatic  environment.  Mining  operations  and  metal  refining  and  processing 
industries  are  important  point  sources  of  trace  metals.  In  certain  parts  of  the  world  (Eastern 
Canada) ,  atmospheric  deposition  also  represents  an  important  direct  contribution  of  trace 
metals  to  surface  waters;  the  concomitant  acidic  deposition  may  also  lead  to  increased 
mobilization  of  trace  metals  from  terrestrial  sources  to  the  aquatic  system.  Asa.  result  of 
complex  physical,  chemical  and  biological  processes,  most  of  the  trace  metals  introduced  into 
the  aquatic  environment  are  found  associated  with  the  bottom  sediments,  where  they  constitute 
a  potential  danger  for  benthic  organisms. 

Management  of  in-place  contaminants  is  one  of  the  major  problems  confronted  by 
governmental  agencies  responsible  for  the  protection  of  the  environment.  In  principle,  remedial 
actions  such  as  reduction  in  waste  disposal,  chemical  treatment  of  in-place  contaminants, 
capping  or  dredging  can  be  undertaken;  these  actions  are  however  costiy.  It  is  clear  that  the 
development  of  rational,  effective  and  economical  strategies  to  solve  the  problem  posed  by 
sediment-bound  toxic  metals  will  depend  greatiy  on  our  ability  to  predict  how  remedial  actions 
will  improve  water  quality  and  how  these  changed  conditions  will  affect  aquatic  organisms. 
Important  progress  in  this  direction  wiU  be  made  only  by  understanding  the  biogeochemical 
processes  governing  metal  accumulation  by  benthic  organisms  under  field  situations. 


The  general  objective  of  the  present  project  was  to  develop  and  test  in  the  field  a  deterministic 
model,  based  on  well  accepted  theories,  to  predict  metal  accumulation  in  benthic  organisms. 
The  idea  behind  this  development  was  that  such  a  model,  because  of  its  mechanistic  basis, 
should  be  of  general  applicability  (i.e.  it  should  apply  to  sites  others  than  those  used  in  the 
model  calibration).  The  specific  objectives  were:  i)  to  derive  predictive  equations  for  trace  metal 
concentrations  in  bivalves  that  are  valid  for  any  aquatic  environment  where  these  organisms 
can  be  found;  ii)  to  evaluate  the  utility  of  the  bivalves  as  bioindicators  for  trace  metals.  By 
coupling  concepts  from  the  free-ion  model  of  trace  metal  organisms  interactions  with  the 
surface  complexation  theory,  we  have  made  substantial  progress  towards  this  general 
objective,  using  Cd  as  our  test  metal. 

2.  MATERIAL  AND  METHODS. 

2.1.  Study  area. 

The  sediments,  associated  porewater  samples  and  bivalve  specimens  were  collected  at 
littoral  stations  in  lakes  located  in  the  Muskoka  area,  Ontario  (Table  1;  Figure  1).  They  lie  on  the 
Precambrian  Shield.  The  closest  areas  of  industrial  activities  are  Sudbury  («  2(X)  km  northwest) 
and  Toronto  (~  250  km  south). 

2.2.  Sampling. 

The  interstitial  and  overlying  water  samples  were  obtained  with  in  situ  samplers  (porewater 
peepers;  1  cm  vertical  resolution;  3  per  lake)  similar  to  those  described  by  Hesslein  (1976)  and 
Carignan  et  al.  (1985).  The  peepers  were  filled  with  demineralized  water  and  bubbled  with 
nitrogen  for  at  least  24  h  in  plexiglass  cylinders  filled  with  demineralized  water  before  being 
inserted  vertically  in  the  lake  sediments  by  a  diver.  After  a  two  week  equilibration  period,  the 
peepers  were  retrieved  from  the  sediments  by  a  diver;  samples  were  collected  immediately 
from  five  compartments  above  and  five  below  the  water-sediment  interface,  and  pH  was 
measured  in  each  sample  (Carignan,  1984).  Samples  for  dissolved  sulfate  and  chloride  analysis 
were  removed  from  the  peepers  with  a  syringe  and  injected  into  prewashed  polypropylene 
tubes;  those  for  inorganic  carbon  determinations  were  obtained  with  a  syringe  and  injected 
through  a  septum  into  pre-evacuated  and  pre-washed  glass  tubes.  The  samples  for  the 
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Figure  1 .  Location  of  sampling  stations  in  the  Muskoka  area. 


Table  1.  Sample  identification  and  location  of  sampling  stations  in 
the  Muskoka  area. 

Site       Lake  Location 

I.  D. 


BI-01 

Bird 

45°02' 

'N 

79004 

'W 

PA-01 

Paint 

45»13 

'N 

78°57 

'W 

RI-01 

Ril 

45<'10 

'N 

79000 

•w 

TO- 01 

Tock 

45«16 

'N 

78''53 

•w 

WR-01 

Little 

Wren 

45°11 

'N 

78°51 

•w 

analysis  of  trace  metals  were  collected  by  piercing  the  membrane  (Gelman  HT-200,  0.2  ij.m 
nominal  pore  size)  with  a  Gilson  pipette  fitted  with  acid  cleaned  tips:  these  samples  were 
injected  into  pre-washed  and  pre-acidified  (30  uL  of  IN  Ultrex  HNO3,  final  pH  <  2.5)  Teflon'* 
vials.  Additional  water  was  collected  from  the  compartments  6  and  7  cm  above  the  sediment- 
water  interface  of  each  peeper  (combined  content  5  mL)  for  analysis,  if  necessary,  by  multiple 
injection  in  flameless  atomic  absorption  spectrophotometry  (see  below).  On  several  occasions 
during  the  summer,  water  samples  were  also  collected  by  divers  at  many  of  the  sites  with  a 
clean  polyethylene  bottle  close  to  the  sediment-water  interface  for  pH  measurement.  The 
purpose  of  this  measurement  was  to  calculate  time-averaged  pH  values  (-log  (sCH'^J/n))  for  the 
summer  period. 

Sediment  cores  (4-5  per  site)  were  collected  by  a  diver,  close  to  the  peepers,  with  plexiglass 
tubes  (9  cm  diameter).  The  tubes  were  tightly  closed  to  minimize  perturbations  of  the  sediments 
during  their  transport  to  the  shore.  The  sediment  cores  were  extruded  on  the  shore,  and  only 
the  uppermost  half  centimeter,  containing  oxidized  sediments,  was  retained.  These  samples 
were  placed  in  500  mL  centrifugation  bottles  half  filled  with  lake  water,  and  kept  at  !>i4°C  during 
transport  to  the  laboratory  where  they  were  kept  frozen  until  analysis. 

Specimens  oiElliptio  complanata  (10  per  site)  were  obtained  by  a  diver  at  each  site,  within  a 
radius  of  =  50  m  from  the  sediment  and  porewater  collection  site.  The  bivalves  were  placed  in 
plastic  bags  with  lakewater  and  maintained  at  ~  4  °C  during  transport  to  the  laboratory  where 
they  were  left  to  depurate  for  at  least  24  h  in  aerated  lake  water.  They  were  then  dissected  into 
gills,  mantle,  hepatopancreas  and  remaining  tissues,  hereafter  referred  to  as  "remains".  For 
each  sampling  site,  individual  tissue  types  from  the  ten  animals  collected  were  pooled  and 
frozen  at  -20  °C  until  needed  for  analysis. 

2.3  Analyses. 

The  water  samples  were  analyzed  for  Cd,  Cu,  Ni,  Pb,  Zn,  Fe,  Mn,  SO^,  CI,  major  cations  and 
organic  and  inorganic  carbon.  The  metal  concentrations  were  obtained  by  flame  atomic 
absorption  spectrophotometry  (AA)  when  possible  (Fe,  Mn,  Ca,  Mg,  Na,  K;  Varian  Techtron, 
Model  575ABQ  or  Model  Spectra.  AA-20)  or,  otherwise,  by  flameless  A  A  (Fe,  Mn,  Cd,  Cu,  Ni, 
Pb,  Zn;  Varian  Techtron,  Model  1275  or  Spectra  AA-30;  GTA-95  or  GTA-96).  When  initially 


measured  levels  of  cadmium  in  the  overlying  waters  did  not  exceed  a  concentration  threshold 
of  five  times  the  analytical  detection  limit  (0. 13  nmol  L"'),  its  concentration  was  determined  by  a 
multiple  injection  technique.  The  necessary  number  of  aliquots  were  injected  successively  and 
subjected  to  the  preliminary  drying  steps;  the  combined  sample  was  then  atomized.  The 
detection  limits  for  Ni  and  Pb  were  not  improved  with  this  technique.  The-National  Research 
Council  of  Canada  (NRC)  Riverine  Water  Reference  Material  (SLRS-2)  was  analysed  for  Cd,  Cu 
and  Zn;  we  obtained  0.024±0.(X)3  ng  Cd  L'^  (N=4;  certified  value  0.028±0.004  Mg  Cd  L'), 
2.96±0.13  Mg  Cu  L-'  (N=6;  certified  value  2.76±0.17  Mg  L"*)  and  3.38±0.07  Mg  Zn  L'  (N=7; 
certified  value  3.33±0.15  ng  L').  The  National  Bureau  of  Standards  Reference  Material  1643b 
(Trace  Elements  in  Water)  was  analyzed  after  lOX  to  40X  dilutions  for  Pb;  we  obtained: 
25.7±2.4  Mg  Pb  L-'  (N=4;  certified  value  24.1  ±0.7  ixg  Pb  L"').  Sulfate  and  chloride 
concentrations  were  determined  by  ion  chromatography  (Dionex  Autolon,  System  12); 
dissolved  inorganic  carbon  was  measured  by  gas  chromatography  (Carignan  1984)  and 
dissolved  organic  carbon  by  persulfate-ultraviolet  oxidation,  followed  by  conductometric 
determination  of  the  released  COj  on  a  Technicon  autoanalyser. 

The  surficial  sediment  samples  were  thawed  and  centrifuged  to  remove  excess  water; 
subsamples  (equivalent  to  a  1  g  dry  wt.)  were  extracted  and  trace  metals  were  partitioned  into 
the  following  empirical  fractions  (Tessier  et  al.  1985):  i)  the  sediment  subsample  was  extracted 
for  30  min  with  1  N  MgCl2;  ii)  the  residue  from  (i)  was  extracted  for  5  h.  with  CHjCOONa 
adjusted  to  pH  5  with  CH3COOH;  iii)  the  residue  from  (ii)  was  extracted  for  30  min.  at  room 
temperature  with  0. 1  M  NHjOH  •  HCl  in  0. 1  N  HNO3;  iv)  the  residue  from  (iii)  was  extracted  for  6 
h.  at  96°C  with  0.04  M  NH^OH-HCl  in  25%  (v/v)  CH3COOH;  v)  the  residue  from  (iv)  was 
extracted  for  5  h.  at  85° C  with  B^O^  adjusted  to  pH  2  with  HNO3,  and  then  at  room  temperature 
with  3.2  M  NH^OAc  in  20%  (v/v)  HNO3;  vi)  the  residue  from  (v)  was  digested  successively  with 
concentrated  HNO3  (15  min;  under  reflux;  followed  by  evaporation  to  dryness),  concentrated 
HCLO4  ^^  ^^''  ^  ^^^  under  reflux)  and  HF  (15  mL;  evaporation  to  dryness;  12-18  h)  and 
dissolution  of  the  residue  was  effected  in  50  mL  of  5%  HCL.  Fe,  Mn,  Cd,  Cu,  Ni,  Pb  and  Zn 
concentrations  in  the  extracts  were  determined  by  flame  or  flameless  AAS,  using  the 
appropriate  extractant  matrices  for  standards  and  blanks.  Details  of  these  procedures  are  given 
in  earlier  publications  (Tessier  et  al.  1979,  1985).  The  National  Research  Council  of  Canada 
reference  sample  of  marine  sediment  MESS-1  has  been  subjected  to  the  sequential  extraction 


procedure.  We  measured,  for  the  sum  of  the  six  fractions:  25.2+1.4  ng  Cu  g'  (N=4;  certified 
value  25.1  ±3.8  ng  Cu  g'),  169±2  Mg  Zn  g"'  (N=4;  certified  value  191  ±17  ng  Zn  g"'), 
29.6±0.9  Mg  Pb  g'  (N=4;  certified  value  34.0±6.1  Mg  Pb  g"').  For  Cd,  and  Ni,  which  showed 
lower  than  detection  limit  values  in  some  fractions,  we  found  the  following  ranges 
0.56±0.07  <  X  <  0.72±0.03  Mg  Cd  g'  (N=4;  certified  value  0.59±0.10  Mg  Cd  g"')  and 
26.0±1.5  <  X  <  28.6±2.8  Mg  Ni  g*'  (N=4;  certified  value  29.5±2.7  Mg  Ni  g').  Sediment  organic 
carbon  concentrations,  {C    },  were  determined  with  a  CNS  analyzer  (Carlo-Erba,  Model 
NA1500)  after  removal  of  inorganic  carbon  by  acidification  with  diluted  H2SO4  (0.5  mol  L"';  15 
min.,  100  mL/g  sediment  dry  weight). 

Each  pooled  bivalve  tissue  was  homogenized  (Brinkman  tissue  grinder,  Model  CH-6010).  A 
subsample  was  then  dried  to  constant  weight  to  determine  the  wet:dry  weight  ratio,  and  a 
second  subsample  was  digested  in  a  Teflon  bomb  with  concentrated  nitric  acid  (Aristar; 
3  mL  g"'  of  tissue  dry  weight)  in  a  microwave  oven  at  pressures  between  5500  and  7000  KPa  for 
a  1  min.  A  certified  reference  material  (lobster  hepatopancreas,  TORT-1,  NRC)  was  submitted 
to  the  same  digestion  procedure  and  analyzed  for  Cd,  Cu,  Ni,  Pb  and  Zn;  we  obtained 
24.1  ±0.5  Mg  Cd  g-'  (N=5;  certified  value  26.3±2.1  Mg  Cd  g'*),  385±5  Mg  Cu  g"'  (N=5;  certified 
value  439±22  Mg  Cu  g''),  2.1  ±0.3  Mg  Ni  g'^  (N=4;  certified  value  2.3±0.3  Mg  Ni  g''), 
9.6±0.7  Mg  Pb  g'  (N=5;  certified  value  10.4±2.0  Mg  Pb  g-*)  and  182±31  Mg  Zn  g'  (N=5; 
certified  value  177±10  Mg  Zn  g'). 


3.  RESULTS 

The  concentrations  of  dissolved  constituents  above  the  sediment-water  interface,  as 
obtained  with  the  porewater  peepers,  are  given  in  Table  2.  Dissolved  cadmium  concentrations 
vary  between  0.04  and  0.25  nmol  L';  these  values  are  much  lower  than  those  reported  by 
Hinch  and  Stephenson  (1987)  for  Beech  (3  nmol  L"')  and  Tock  (3  nmol  L')  Lakes.  For  this  latter 
lake,  we  measured  0.23  nmol  L"'.  Campbell  and  Evans  (1991)  measured  values  between  0.2 
and  1.25  nmol  L"'  in  the  epilimnia  of  various  lakes  in  south-central  Ontario;  some  of  their 
sampled  lakes  are  in  the  Muskoka  area.  We  measured  dissolved  Cu  concentrations  between 
4.6  and  7.0  nmol  L''  and  Zn  concentrations  between  17  and  89  nmol  L"';  these  values  are 


Table  2.  Mean  concentrations  (±SD)  of  total  dissolved  constituents 
(Ca,  Cd,  Cl,  Cu,  Fe,  Mn,  Ni,  Pb,  SO4 ,  Zn,  pH,  inorganic  and  organic 
carbon)  above  the  sediment-water  interface,  as  obtained  with  the 
porewater  peepers. 


Site 

PH 

Cd 

Cu 

Ni 

Pb 

Zn 

( 

nmo 

1  L--L) 

BI-01 

6.37 

0.17^ 

5.5±1. 

,1 

<6 

<1 

44.4±5.9 

PA- 01 

7.13 

0.04310.006 

7.0±0. 

.8 

<6 

<1 

16-7±2.4 

RI-01 

6.02 

0.20±0.03 

5.2±0. 

.8 

<6 

<1 

89.2±5.1 

TO- 01 

6.43 

0.23^ 

6.2^. 

<6 

<1 

71.2113.7 

WR-01 

6.27 

0.25±0.01 

4.6±1, 

.1 

<6 

<1 

80.8110.1 

Site 

Ca 

Mg 

Na 

K 

CI 

SO4 

(Mmol  L"-'-) 

BI-01 

107137 

34.611, 

.2 

3517 

1213 

1511 

7816 

PA- 01 

147136 

51.410, 

.6 

8615 

1612 

8313 

9713 

RI-01 

113111 

35.012, 

.1 

4410 

1010 

3512 

7814 

TO-01 

139140 

30.811, 

.2 

3910 

1010 

4113 

7913 

WR-01 

80117 

28.114, 

.6 

6417 

1413 

5816 

8117 

Site 

Fe 

Mn 

^inorg 

^org 

(Minol  L"l) 

(/imol 

C  L~l) 

BI-01 

1.310.9 

0, 

.4510. 

,13 

210160 

4751118 

PA- 01 

0.910.8 

0, 

.1110. 

,08 

130130 

440134 

RI-01 

1.810.2 

0, 

.7510. 

,15 

180120 

473160 

TO-01 

0.510.3 

0, 

.3210. 

.11 

110150 

474199 

WR-01 

1.310.8 

0, 

.6010. 

,14 

110180 

453163 

mean  of  two  values. 


slightly  lower  than  those  reported  by  Hinch  and  Stephenson  (1987)  for  Beech  (16  nmol  Cu  L'; 
77  nmol  Zn  L')  and  Tock  (16  nmol  Cu  L'';  120  nmol  Zn  L')  Lakes.  The  concentrations  of 
dissolved  Ni  and  Pb  were  below  our  detection  limits  (  6  and  1  nmol  L'  respectively).  For  these 
last  two  metals,  multiple  injections  did  not  improve  detection  limits. 

Concentrations  of  sedimentary  trace  metals  obtained  in  various  fractions  vary  slightly  from 
one  lake  to  another  (Table  3).  Comparison  of  sediment  trace  metal  concentrations  among  lakes 
or  with  those  obtained  by  other  researchers  is  difficult,  since  the  sediment  samples  were 
obtained  from  the  littoral  zone  of  the  lakes,  where  sediment  composition  is  often  patchy.  The 
results  show,  however,  general  trends  that  were  observed  in  other  studies  (Forstner  and 
Wittmann,  1981).  For  example,  non-residual  Cu  is  found  mostly  in  fraction  F4  (between  65% 
(Paint  Lake)  and  97%  (Little  Wren  Lake)),  which  confirms  the  high  affinity  of  natural  organic 
matter  for  Cu.  Large  proportions  of  non-residual  Pb  (between  47%  (Tock  Lake)  and  87%  (Paint 
Lake))  and  non-residual  Zn  (between  46%  (Ril  Lake)  and  64%  (Paint  Lake))  are  found  in  the 
sum  of  fractions  F3A  and  F3B;  similar  results  were  reported  for  suspended  river  sediments 
(Tessier  et  al,  1980).  The  ratio  {Cj,^}:{Fe-ox}  in  the  sediments  ({Fe-ox}  is  the  sum  of 
{Fe(F3A)}  and  {Fe(F3B)}  in  Table  3)  varies  between  10  (Paint  Lake)  and  130  (Little  Wren 
Lake);  this  large  variation  in  the  ratio  explains  at  least  partly  variations  in  metal  partitioning 
observed  among  lakes.  For  example,  the  ratio  [Cu(F4)]:([Cu(F3A)]  -I-  [Cu(F3B)])  is  linearly 
correlated  (r-=0.92;  N=5)  with  the  ratio  {C}: {Fe-ox}.  This  strong  correlation  suggests  a 
competition  between  iron  oxyhydroxides  and  organic  matter  for  copper;  such  a  competition 
has  been  reported  by  Luoma  (1986). 

Concentrations  of  Cd,  Cu,  Ni,  Pb  and  Zn  in  the  tissue  of  £.  complanata  are  given  in  Table  4.  We 
found  that  cadmium  concentrations  in  the  tissues  of  the  bivalve  varied  between  13  and 
26  Mg  g"'-  These  values  are  similar  to  those  reported  for  E.  complanata  by  Servos  et  al. ,  (1987) 
for  Lake  of  Bays  (12  Mg  g'')  and  by  Hinch  and  Stephenson  (1987)  for  Beech  (1 1  Mg  g')  and 
Tock  (15  Mg  g')  Lakes.  The  values  obtained  for  Cu  (between  5.9  and  7.3  Mg  g"')  and  Zn 
(between  170  and  215  Mg  g"')  are  similar  to  those  reported  by  Hinch  and  Stephenson  (1987)  for 
Beech  (7.3  Mg  Cu  g"';  155  Mg  Zn  g"')  and  Tock  (10  Mg  Cu  g"';  130  Mg  Zn  g"')  Lakes. 
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Table  3.  Partitioning^*)  of  trace  metals  in  the  oxic  layer  of  lake 
sediments  and  concentration  of  sedimentary  organic  carbon.  The 
concentrations  are  given  on  a  dry  weight  basis. 

Site     {M(F1)}   {M(F2)}   {M(F3)}   {M(F3B) }   {M(F4) }   {M(F5)}   {M}^ 

{Cd}  (nmol  g~^) 


BI-01 

1.12 

1.62 

0.89 

4.80 

<0.5 

<0.9 

8.43 

PA- 01 

0.61 

0.82 

0.36 

0.18 

0.71 

<0.9 

2.67 

RI-01 

1.00 

1.77 

1.33 

5.16 

0.53 

<0.9 

9.79 

TO-01 

0.90 

3.61 

1.69 

10.5 

1.25 

<0.9 

17.9 

WR-01 

0.71 

1.05 

0.80 

7.38 

1.25 

<0.9 

11.2 

{Cu}  (nmol  g"^) 


BI-01 

1.97 

0.27 

1.42 

5.82 

79.6 

49. 

,6 

139 

PA- 01 

<0.5 

2.77 

2.05 

8.50 

24.2 

26. 

.8 

64. 

3 

RI-01 

0.87 

1.87 

1.42 

5.35 

79.9 

49. 

,6 

139 

TO-01 

<0.5 

5.52 

1.26 

11.3 

187 

125 

331 

WR-01 

<0.5 

<0.2 

<0.2 

2.68 

94.7 

27. 

.5 

125 

{Ni}  (nmol  g  ^) 


BI-01 

6.88 

7.63 

<7 

35.4 

51.1 

162 

263 

PA- 01 

3.99 

7.63 

9 

25.9 

10.2 

145 

201 

RI-01 

21.3 

6.82 

<7 

32.7 

34.1 

153 

248 

TO-01 

23.2 

13.1 

<7 

72.4 

129 

162 

400 

WR-01 

6.93 

<7 

2.7 

39.0 

64.7 

145 

258 

{Pb}  (nmol  g"^) 


BI-01 

1.88 

27.5 

2.08 

112 

18. 

.3 

57. 

.9 

220 

PA- 01 

<1 

2.7 

4.05 

33, 

,1 

2. 

.90 

55, 

.5 

98. 

3 

RI-01 

16.7 

34.2 

2.08 

85. 

,9 

14, 

.5 

57, 

.9 

211 

TO-01 

6.56 

79.1 

2.99 

138 

75, 

.3 

53, 

.1 

355 

WR-01 

3.38 

19.8 

<2 

48. 

,7 

26, 

.1 

45, 

.8 

144 

{Zn}  (nmol  g  ^) 


BI-01 

177 

173 

78.5 

503 

97.9 

667 

1697 

PA- 01 

43.7 

108 

69.7 

283 

47.7 

531 

1083 

RI-01 

358 

182 

124 

404 

82.6 

541 

1691 

TO-01 

203 

378 

190 

1361 

283 

765 

3180 

WR-01 

110 

145 

82.7 

491 

135 

505 

1469 

Table    3.     (Continue). 

Site  {M(F1)}       {M(F2)}       {M(F3)}       {M(F3B)}       {M(F4)}       {M(F5)}       {M}^ 


{Fe}  (/imol 

g" 

■h 

BI-01 

0. 

.20 

8, 

.49 

6. 

.59 

95.8 

10. 

.3 

413 

534 

PA-01 

0. 

.20 

4, 

.03 

6, 

.52 

105 

3. 

.65 

397 

517 

RI-01 

0. 

.04 

8, 

.83 

9, 

.79 

73.9 

10. 

.6 

380 

483 

TO-01 

0. 

.06 

6, 

.79 

7. 

.07 

113 

33. 

.9 

544 

705 

WR-01 

0. 

,14 

7, 

.72 

7. 

.95 

63.8 

24. 

.9 

344 

448 

{Mn}    (/Ltmol  g~^) 


BI-Ol 

0.87 

1.26 

2.29 

1.49 

0.06 

9.56 

15.5 

PA-01 

0.55 

2.02 

25.3 

3.08 

0.06 

9.65 

40.7 

RI-01 

0.69 

1.24 

1.97 

1.04 

0.05 

7.88 

12.9 

TO-01 

0.88 

1.48 

0.87 

1.18 

0.05 

6.19 

10.7 

WR-01 

1.07 

0.35 

0.07 

0.28 

0.08 

9.08 

10.9 

^^org}.-. 
(Mitiol   C   g    ^) 


BI-Ol  4472 

PA-01  1107 

RI-01  3230 

TO-01  9967 

WR-01  9358 


{M (FI)}  {M(F5) }  represent  metal  concentrations  in  fractions  1 

to  6  following  the  sequence  given  in  the  text. 


Table  4.  Metal  concentrations  (on  a  dry  weight  basis)  in  the  tissues  of  Elliptio  complanata  from 
various  sampling  sites.  He,  hepatopancreas;  Gi,  gills;  Ma,  mantle;  Re,  remaining;  org,  whole 
organisms. 


Site 

[M(He)] 

[M(Gi 

•)] 

Cd 

[M(Ma)] 

(Mg  g"^) 

[M(Re)] 

[M(org)] 

Bl-dl 

9.40 

7.54 

16.8 

12.9 

12.7 

PA- 01 

17.9 

8.28 

14.9 

7.59 

10.3 

RI-01 

17.3 

14.0 

24.6 

12.3 

15.4 

TO-01 

25.9 

18.1 

47.8 

20.3 

25.7 

WR-01 

14.6 

17.1 

31.6 

10.0 

15.7 

Cu    (/ig  g"^) 


BI-01 

10.7 

13.2 

5.86 

5.37 

7.14 

PA- 01 

12.2 

11.8 

6.58 

5.31 

7.29 

RI-01 

10.5 

13.4 

4.92 

4.82 

6.63 

TO-01 

12.6 

12.8 

6.35 

4.75 

7.32 

WR-01 

8.05 

10.2 

3.68 

Ni  (Mg  g"^) 

4.85 

5.89 

BI-01 

1.59 

2.59 

1.46 

1.09 

1.40 

PA- 01 

0.82 

2.88 

2.14 

0.77 

1.34 

RI-01 

1.68 

4.81 

3.31 

1.82 

2.45 

TO-01 

2.47 

2.89 

2.46 

1.38 

1.94 

WR-01 

1.39 

2.86 

2.75 

pb  (Mg  g"^) 

0.92 

1.63 

BI-01 

42.0 

14.0 

51.1 

37.2 

38.8 

PA- 01 

1.81 

1.26 

1.97 

1.54 

1.61 

RI-01 

20.7 

10.5 

21.7 

13.0 

15.3 

TO-01 

92.4 

37.4 

102 

53.4 

65.3 

WR-Ol 

27.8 

12.2 

65.7 

zn  (Mg  g"-'-) 

25.0 

30.7 

BI-01 

168 

234 

197 

153 

173 

PA- 01 

128 

289 

266 

120 

173 

RI-01 

158 

356 

254 

135 

187 

TO-01 

194 

330 

309 

159 

215 

WR-01 

144 

282 

257 

120 

175 

4.  DISCUSSION 

4.1  The  free-ion  activity  model. 

According  to  the  free-ion  activity  model  (Figure  2),  the  effects  (including  bioaccumulation)  of 
a  trace  metal  M  on  an  organism  which  obtains  the  trace  metal  from  the  water  should  be 
predicted  by  the  free  ion  concentration  (or  activity)  of  this  metal,  [M^"^].  The  concentration  of  a 
metal  accumulated  into  an  organism,  [M(org)],  should  thus  be  written  (Morel,  1983): 

[M(org)]      =      f([MZ+])  ^  (1) 

where  f  means  a  function.  This  model  has  been  developed  for  explaining  (with  great  success) 
biological  effects  of  trace  metals  on  unicellular  organisms.  Laboratory  bioassays,  conducted 
with  benthic  animals  exposed  to  dissolved  metals  under  carefully  controlled  conditions,  have 
shown  that  the  effects  of  trace  metals  on  these  animals  are  also  related  to  the  free  metal  ion 
concentration  and  not  to  the  total  metal  concentrations  in  the  animals'  environment.  For 
example,  mortality  of  the  shrimp  Palaemonetes  pugio  exposed  to. cadmium  was  shown  to  be 
related  to  [Cd^"^]  (Sunda  et  al,  1978).  Similarly,  short-term  (14  d)  accumulation  of  copper  in  the 
American  oyster  Crassostrea  virginica  was  found  to  be  related  to  [Cu^"^]  (Zamuda  and  Sunda, 
1982).   Note,  however,  that  the  applicability  of  the  free-ion  activity  model  in  the  field,  where 
dissolved  ligands  (e.g.  natural  organic  matter)  differ  from  the  synthetic  ligands  used  in  the 
laboratory  (e.g.  NT  A,  EDTA),  has  never  been  demonstrated. 

For  the  following  discussion,  we  have  focused  on  cadmium;  among  the  metals  measured, 
water,  sediment  and  animal  data  are  more  complete  for  this  metal  (other  metals  show  many 
"lower  than  detection  limit"  values,  especially  for  the  water  samples).  To  enhance  the 
representativity  of  the  study,  we  include  also  Cd  data  obtained  for  other  lakes  that  we  have 
studied,  which  are  located  in  the  areas  of  Sudbury  (On),  Rouyn-Noranda  (Qc;  mining  area; 
heavily  contaminated  with  Cd),  Chibougamau  (Qc;  mining  area),  Québec  (Qc)  and  Eastern 
Townships  (Qc).  These  lakes  are  distributed  over  a  geographical  area  of  about  350  000  km^ 
and  they  represent  a  large  gradient  of  Cd  contamination  and  a  wide  range  of  pH  values.  Since 
additional  measurements  are  needed  to  obtain  a  satisfactory  picture  of  Cd  concentrations  in  E. 
complaruua  (see  below;  this  section  and  section  4.3),  we  also  include  in  the  following 
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Figure  2.  Schematic  representation  of  the  various  reactions  occurring  in  the  water  overlying  the 
sediments.  M'*  is  the  free  metal  ion;  L^  and  L^^  are  dissolved  inorganic  and  organic  ligands 
respectively;  L  is  a  ligand  at  the  surface  of  an  organism  or  a  transport  molecule;  S  is  a 
component  of  tde  sediment  (e.g.  iron  oxyhydroxides;  organic  matter)  which  sorbs  tne  trace 
metal;  ML^^^,  ML^^^,  ML^  and  MS^,  represent  the  metal  bound  to  the  various  ligands  or  sediment 
components. 
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discussion  Cd  data  that  we  have  obtained  for  another  widespread  freshwater  bivalve,  Anodonta 
grandis. 

We  have  plotted  in  Figure  3  cadmium  concentrations  in  E.  complanata  (Fig.  3a)  and  in  A. 
grandis  (Fig.  3b)  as  a  function  of  [Cd^^].  We  have  calculated  [Cd^+]  with  the  equilibrium  model 
HYDRAQL  (Papelis  et  al. ,  1988),  using  the  total  dissolved  Cd  concentrations  ([Cd];  Table  2), 
the  measured  concentrations  of  the  inorganic  ligands  concentrations  (Table  2),  and  the  stability 
constants  of  the  inorganic  complexes  (Smith  and  Martell,  1977);  the  stability  constants  used  for 
calculating  Cd  complexation  are  given  in  Table  5.  Because  the  thermodynamic  equilibrium 
constants  for  Cd-DOM  (dissolved  organic  matter)  are  unknown,  possible  cadmium 
complexation  by  natural  organic  ligands  could  not  be  considered  in  the  calculation. 

According  to  the  data  in  Figure  3b,  Cd  concentrations  in  the  tissues  of  ^.  grandis  (/xg  g"')  can 
be  described  by  : 

[Cd(org)]      =      44    (±5)     [Cd^+J    +    10    (±15)  (2) 

(n  =   17;    r^   =   o.81;    [Cd^^]    in  nmol   L"^) 

The  relationship  between  cadmium  concentration  in  E.  complanata  and  [Cd^"^]  is,  however,  not 
as  clear;  if  data  from  D'Alembert  Lake  are  removed,  it  would  resemble  an  hyperbolic  function. 
Additional  measurements,  especially  at  sites  where  dissolved  Cd  concentrations  are  high,  are 
needed  before  a  definitive  interpretation  can  be  made. 

It  might  be  argued  that  much  of  the  scatter  in  Figure  3a,b  is  due  to  variation  in  dissolved  Cd 
on  short  time  scales  (Tessier  et  al. ,  1993).  For  example,  Yan  et  al.  (1990)  reported  a  2-fold 
change  in  dissolved  Cd  concentration  in  Red  Chalk  Lake  during  the  ice-free  season;  this  lake  is 
located  in  the  Muskoka  area.  Dissolved  Cd  concentrations  used  in  the  present  study  are  mean 
values  for  a  single  sampling.  Cadmium  concentrations  in  the  bivalves  are  probably  not  sensitive 
to  short-term  variations  in  dissolved  Cd.  For  example,  we  have  observed  that  Cd 
concentrations  in  the  tissues  of  A.  grandis  specimens  transplanted  from  the  unpolluted  Lake 
Brompton  (Eastern  townships)  to  the  polluted  Lake  Joannes  (Rouyn-Noranda  area)  had  not 
reached,  after  three  years,  the  cadmium  levels  observed  in  the  indigenous  specimens  from 
Lake  Joannes.  The  use  of  time-averaged  values  of  dissolved  Cd,  instead  of  the  single  values 


(^.6  6r()  [(6jo)po] 


•  S 

ë  ^ 

<  g 

s  s 

^    • 

~-    en 

Ci    ™ 

«*"  S 
o  a 

eg 

o-o 
oU 


•S3 


o 


■o 
O 


(^.6  6t1)  [(5jo)po] 


.s  a 
•Sl 

c  o 

8§ 

11 
II 


c  o 

l.i 

IS 


;:  E  4) 
J-  ca  «n 

o 

c 


C3£ 


rt    O  -^ 


^  îa  o 

PU   w  > 


Table  5.  Stability  constants  used  in  the  equilibrium  model  HYDRAQL 
(Papelis  et  al.,  1988)  for  the  complexation  with  inorganic  ligands. 

Reaction  Log  K 

Cd^"*"   +  C02^~      =         CdCOj  3.20 

Cd^"^   +   SO42-   =    CdS04  2.46 

Cd2+   +   Cl"     =    CdCl'*'  2.00 

Cd2+   +   2  Cl"   =    CdClj  2.60 

Cd2+   +   3  Cl~   =    CdCl3"  2.40 

Cd2+   +   OH"     =    CdOH"^  -10.10 

Cd^"»"   +   2  OH"   =    Cd(0H)2  -20.40 

Cd2+   +   3  OH"   =    Cd(0H)3"  -33.00 

Cd2+   +   4  OH"   =    Cd(OH)42-  -47.35 


used  here,  would  improve  the  predictive  power  of  Eqn  (2)  and  help  in  choosing  the  proper 
function  to  predict  cadmium  concentration  in  E.  complanata  (Fig.  3a).  Measurement  of  temporal 
changes  in  Cd  concentrations  in  the  tissues  of  transplanted  A.  grandis  specimens  suggests 
that  they  accumulate  metals  only  during  the  warm  season  (e.g.  from  May  to  October).  For 
optimum  predictive  power,  the  mean  [Cd^"*"]  should  thus  be  preferentially  obtained  for  this 
period. 

The  data  shown  in  Figure  3  are  consistent  with  the  free-ion  activity  model.  It  should  be  noted, 
however,  that  very  similar  relations  would  have  been  obtained  with  total  dissolved  Cd  ([Cd]), 
since  its  concentration  at  all  sites  is  close  to  the  computed  value  for  [Cd^"^].  Given  the  similarity 
between  [Cd^"^]  and  [Cd]  values,  and  the  scatter  shown  in  Fig.  3,  it  is  not  possible  with  the 
present  data  to  demonstrate  unambiguously  that  the  free-ion  activity  model  applies  in  the  field 
to  Cd  accumulation  in  E.  complanata  or  A.  grandis.  To  do  so,  one  would  need  a  natural  setting 
where  the  degree  of  complexation,  a  =  [Cd]/[Cd^"^],  vary  from  site  to  site.  Both  SO^  and  CI  are 
known  to  complex  Cd,  but  it  is  unlikely  that  E.  complanata  and  A.  grandis  would  be  found  in 
freshwaters  where  these  inorganic  ligands  are  present  in  sufficiently  high  concentrations  to 
complex  Cd  significantly.  An  alternative  approach  to  demonstrate  the  application  of  the  free-ion 
activity  model  to  Cd  accumulation  by  the  two  bivalves  in  the  field  would  be  to  find  sites  where 
natural  organic  matter  complexes  Cd  significantly,  provided  of  course  that  we  can  account  for 
this  complexation  and  correctly  estimate  [Cd^"^]. 


4.2  Relationship  between  dissolved  and  surficial  sedimentary  Cd. 

Cadmium  concentrations  in  the  surficial  oxic  sediments  are  related  to  dissolved  Cd 
concentrations  in  the  overlying  water,  i.e.  in  the  bivalves'  aqueous  environment.  The  partitioning 
of  Cd  between  water  and  surficial  oxic  sediments  has  been  interpreted,  using  surface 
complexation  concepts,  in  terms  of  the  sorption  of  this  metal  to  the  two  following  sediment 
components:  organic  matter  and  iron  oxyhydroxides.  The  development  of  this  partitioning 
model  is  discussed  in  detail  in  the  paper  by  Tessier  et  al. ,  (1993)  given  in  the  Appendix. 
According  to  this  partitioning  model,  [Cd^"^]  can  be  calculated  from  the  following  three 
expressions: 
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{Fe-Cd}     [H"^]^ 

[Cd2  +  ]       =      (3a) 

Npe-  Kpe-Cd  {Fe-ox} 

{OM-Cd}  [H''"]^ 

[Cd2  +  ]   =   (3b) 

NoM-*KoM-Cd  {OM} 

{Cd}T  [H+]^+y 

[Cd2  +  ]   =   (3c) 

NFe-*KFe-Cd  {Fe-ox}  [H+l^  +  NoM-*KoM-Cd  {O^}  [H+]^ 


where: 

-  (Fe-ox}  is  the  concentration  of  iron  oxyhydroxides  in  the  surficial  sediments  (nmol  g"') 
obtained  by  extracting  the  sediments  with  a  reducing  reagent); 

-  {OM}  is  the  concentration  of  organic  carbon  (nmol  g"*)  in  the  surficial  sediments; 

-  {Cd}^.  is  the  total  Cd  concentration  (nmol  g')  in  the  surficial  sediments; 

-  {Fe-Cd}  and  {OM-Cd}  are  the  concentrations  of  cadmium  (nmol  g')  associated  with  iron 
oxyhydroxides  and  organic  matter  in  the  surficial  sediments,  respectively; 

-  X  and  y  are  the  average  apparent  numbers  of  proton  released  per  Cd^"^  ion  adsorbed  on 
iron  oxyhydroxides  and  organic  matter  respectively  (Honeyman  and  Leckie,  1986); 

-  Np^  is  the  number  of  moles  of  sorption  sites  of  the  iron  oxyhydroxides  per  mole  of  iron 
oxyhydroxides; 

-  Nqj^  is  the  number  of  moles  of  adsorption  sites  of  the  organic  matter  per  mole  of  organic 
carbon; 

-  'Kp^^j  and  'K^j^^^  are  apparent  overall  equilibrium  constants  for  the  sorption  of  Cd  on 
iron  oxyhydroxides  and  organic  matter  respectively. 

It  should  be  noted  that  the  values  of  the  geochemical  constants  x  (0.82),  y  (0.97),  Np^.*Kpç.cd 
(10' ^°)  and  Nq^-'Ko^^j  (10"^"*)  have  been  determined  experimentally  from  field  geochemical 
data  (see  the  Appendix). 


20 

It  should  also  be  noted  that  any  of  the  three  equations  (3a,  b  or  c)  can  be  used  to  estimate 
[Cd^"^],  depending  on  which  of  the  sedimentary  variables  ({Fe-Cd},  {OM-Cd},  {Fe-ox},  {OM}, 
{Cd},.  )  are  available  and  have  been  measured  accurately  (in  addition  to  water  pH).  Some 
practical  examples  are  given  below: 

-  If  {Cd}^,  {Fe-ox}  and  {OM}  are  available,  any  of  the  three  equations  can  be  used  since 
{Fe-Cd}  and  {OM-Cd}  can  be  calculated  with  the  use  of  the  geochemical  constants  (see 
the  Appendix); 

-  Eqn  (3a)  can  be  used  if  the  ratio  {Fe-Cd}/ {Fe-ox}  is  obtained  by  reductive  dissolution  of 
the  diagenetic  iron  oxyhydroxides  (and  their  associated  Cd)  deposited  on  Teflon  collectors 
inserted  in  sediments  (Belzile  et  al.,  1989);  relatively  pure  iron  oxyhydroxides  are  obtained 
with  this  collection  technique,  since  we  get  rid  of  the  complex  sediment  matrix.  However,  it 
has  been  shown  that  accurate  measurements  of  {Fe-Cd}  are  difficult  to  obtain  by 
extraction  of  a  whole  sediment  sample  with  a  reducing  agent  (Tessier  et  al. ,  1993);  use  of 
Eqn  (3a)  with  data  obtained  in  this  way  is  not  recommended. 

-  If  only  {Cd}^  and  {OM}  values  are  available,  Eqn  (3b)  can  be  used  as  a  first 
approximation,  and  with  circonspection,  by  making  the  approximation  {Cd}^  ~  {OM-Cd}; 
indeed,  the  association  of  Cd  in  oxic  lake  sediments  appears  to  be  dominated  by  its 
interactions  with  sedimentary  organic  matter  (see  the  Appendix). 

-  If  {Cd}^.,  {Fe-ox}  and  {OM}  values  are  available,  Eqn  (3c)  can  be  used 

4.3  Relating  cadmium  concentrations  in  the  bivalves  to  sediment  characteristics. 

We  have  plotted  in  Figure  4  cadmium  concentrations  in  E.  complanata  (Fig.  4a)  and  A. 
grandis  (Fig.  4b)  as  a  function  of  [Cd^"^]  calculated  with  Eqn.  (3c).  According  to  the  data  in 
Figure  4b,  Cd  concentrations  in  the  tissues  of  A.  grandis  (/xg  g'')  can  be  described  by: 

[Cd(org)]  =  59  (±7)  [Cd^+j   +   ii  (±18)  (4) 

(n  =  19  ;  r^  =  0.82;  [Cd^+j  in  nmol  L"^) 
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It  should  be  noted  that  Eqn.  (4)  has  the  same  form  as  Eqn  (2),  i.e.  it  is  a  straight  line: 

[Cd(org)]    =      F    [Cd2+]      +      [Cd(org)]°  (5) 

where  F  is  the  slope  and  [Cd(org)]°  is  the  (small)  interœpt  on  the  y-axis.  The  function 
relating  Cd  concentration  in  E.  complanata  tO'[Cd^"^]  calculated  with  Eqn.  (3c)  is  not  so 
evident  (Fig.  4a).  If  data  from  Ril  and  Dufay  Lakes  were  removed,  it  would  resemble  again  a 
hyperbolic  function  (as  in  Fig.  3a).  Additional  measurements  are  needed  before  a  definitive 
interpretation  can  be  made  concerning  E.  complanata. 

Several  explanations  can  be  invoked  to  explain  scatter  in  Figure  4,  and  the  behavior  of  Ril 
and  Dufay  Lakes  as  outliers.  The  geochemical  model  used  to  partition  sedimentary  Cd  is 
necessarily  crude:  i)  it  assumes  that  no  sediment  component  other  than  Fe  oxyhydroxides 
and  organic  matter  can  sorb  Cd;  ii)  it  assumes  that  total  sedimentary  organic  matter  sorbs 
Cd  -  it  is  conceivable  that  only  certain  organic  fractions  are  responsible  for  sorption,  and 
that  the  composition  of  the  sedimentary  organic  matter  varies  from  lake  to  lake;  iii)  values  of 
the  geochemical  parameters  Np^,  Nq^,  *Kp^_^,j  and  *Kq^_^j  are  considered  here  to  be 
constants,  whereas  they  might  be  expected  to  vary  with  the  nature  of  sedimentary  iron 
oxyhydroxides  and  organic  matter  (Luoma  and  Davis,  1983).  It  is  probable  that  the  nature 
of  these  two  sediment  components  present  in  the  top  half  centimeter  of  the  lake  sediments 
varies  among  stations.  Extraction  of  iron  with  NH^OH-HCl  is  only  moderately  selective  and 
will  solubilize  various  iron  oxyhydroxide  forms;  similarly,  total  organic  sedimentary  carbon  is 
a  very  crude  estimator  of  the  organic  matter  active  in  the  sorption  process.  According  to 
laboratory  experiments  in  well-defined  systems,  the  sorption  constants  should  also  vary 
with  the  density  of  adsorption  (Benjamin  and  Leckie,  1981)  and  with  the  concentration  of 
particles  (DiToro  et  ah,  1986),  both  of  which  vary  among  the  studied  sites.  During  collection 
of  the  surficial  sediment,  the  oxic  layer  may  have  been  contaminated  by  reduced 
sediments,  leading  to  incorrect  estimation  of  {Cd}.^'  {Fe-ox}  and  {OM}  present  in  the  oxic 
surficial  sediment.  Another  source  of  scatter  is  the  variability  in  cadmium  concentration  in 
the  tissues  of  the  bivalves  due  to  biological  factors  (length  or  age,  reproductive  condition). 
Recommendations  are  given  below  (section  5)  to  improve  prediction  of  Cd  (and  other  trace 
metals)  concentrations  in  benthic  organisms. 
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Cadmium  concentrations  in  the  tissues  of  the  freshwater  bivalves  are  much  better  related 
to  the  parameters  appearing  in  the  right-hand  side  of  Eqns  (3)  than  to  total  sedimentary 
cadmium,  {Cd}^  (compare  Fig.  4  and  5).  Indeed,  when  cadmium  concentrations  in  the 
bivalves  are  compared  with  {Cd}^,  no  significant  relationships  emerge  (Fig.  5).  The  same 
conclusion  is  reached  when  Cd  concentrations  in  the  bivalves  are  compared  with  cadmium 
concentrations  found  in  any  sediment  extract,  whether  or  not  the  extractable  Cd  is 
normalized  with  respect  to  sediment  iron  oxyhydroxide  or  organic  carbon  concentrations. 
Such  empirical  normalizations  have  been  advocated  recently  for  the  prediction  of  As,  Cu, 
Hg  and  Pb  bioaccumulation  in  bivalves  (Table  6:  Langston  1980,  1982;  Luoma  and  Bryan, 
1978;  Tessier  et  al. ,  1983;  1984),  but  they  do  not  seem  to  be  adequate  for  predicting  Cd 
accumulation  in  E.  complanata  and  A.  grandis  for  a  large  variety  of  lakes  as  here. 

Prediction  of  [Cd(org)]  with  the  sedimentary  parameters  appearing  in  Eqns.  (3)  is  as 
good  as  the  prediction  with  the  dissolved  [Cd]  and  the  inorganic  ligand  concentrations 
(compare  Figs.  3  and  4).  However,  use  of  sedimentary  measurements  offers  considerable 
advantages.  The  measurement  of  low  dissolved  Cd  concentrations  in  freshwaters  is  difficult 
and  requires  the  use  of  trace  metal-free  techniques;  possibilities  of  inadvertent 
contamination  are  rampant.  An  additional  problem  is  the  temporal  variation  of  dissolved  Cd; 
considerable  effort  would  be  needed  to  obtain  representative  mean  values  of  total  dissolved 
Cd  at  each  sampling  site.  By  comparison,  the  measurement  of  {Cd}^,  {Fe-ox},  {OM}  and 
lake  pH  is  much  simpler. 

An  important  point  to  note  is  that  the  geochemical  constants  x,  y,  Np^  •  "Kp^.^.^  and 
^OM**^M-cd  determined  in  this  study  are  specific  to  cadmium.  These  types  of  constants 
should  be  metal-specific,  and  not  organism-specific.  Once  they  are  determined  accurately 
for  a  given  metal,  they  could  in  principle  be  used  for  predicting  trace  metals  in  any  organism 
that  obtains  its  metal  burden  from  the  ambient  water.  For  a  specific  organism,  the  only 
requirement  would  then  be  to  determine  F  and  [M(org)]°  (as  in  Eqn  5)  from  a  proper 
calibration.  The  geochemical  constants  are  in  fact  used  only  to  estimate  [M^"^],  which  is  the 
same  for  all  the  organisms  at  a  given  site. 
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4.4  Reconciliation  with  empirical  normalization. 

In  earlier  studies,  field  measurements  along  trace  metal  gradients  have  shown  that  the 
prediction  of  metal  concentrations  in  estuarine  and  freshwater  bivalves,  [M(org)],  was  greatly 
improved  when  the  trace  metal  concentration  extracted  from  the  sediments  was  normalized 
with  respect  to  the  iron  oxide  or  organic  matter  content  of  the  sediments  (Table  5).  Hence,  the 
ratios  {Pb}/{Fe}  and  {As)/{Fe},  all  extracted  with  IN  HCl,  and  the  ratio  {Hg  extracted  with 
HN03}/{0M}  were  found  to  be  the  best  predictors  of  lead,  arsenic  and  mercury 
concentrations  respectively  in  the  estuarine  bivalve  Scrobicularia  plana  (Langston  1980,  1982; 
Luoma  and  Bryan  1978).  Similarly,  the  best  predictors  of  Cu  and  Pb  in  the  tissues  of  the 
fireshwater  bivalves  A.  grandis  and  E.  complanata  were  found  to  be  the  ratios  (Cu}/{Fe}  and 
{Pb}/(Fe}  extracted  with  NHjOH-HCl,  a  reducing  reagent  (Tessier  et  al  1983,  1984).  In  other 
words,  according  to  these  studies,  trace  metal  concentrations  in  the  organisms  are  best 
expressed  as: 

{Sn-M} 
[M(org)]    =k +      tM(org)]°  (6) 

where  k  is  a  proportionality  constant;  {S„}  is  the  concentration  of  a  sediment  component  n 
(which  can  be  Fe  oxyhydroxides  or  organic  matter);  {S^-M}  is  the  concentration  of  the  trace 
metal  M  associated  with  sediment  component  n;  [M(org)]°  is  the  intercept  on  the  y-axis  (which 
is  usually  small). 

The  present  study  offers  a  likely  explanation  for  these  findings.  The  studies  of  Langston 
(1980,  1982)  and  Luoma  and  Bryan  (1978)  were  performed  in  estuaries,  where  the  pH  is 
relatively  constant;  the  freshwater  studies  (Tessier  et  al.  1983,  1984)  were  carried  out  in  three 
lakes  located  in  a  restricted  geographical  area  where  the  lake  pH  is  relatively  constant.  In  such 
cases,  [H"^]"  and  [^A^y  become  approximately  constant,  and  combination  of  either  Eqn.  (3a)  or 
(3b)  with  Eqn.  (5)  leads  to  Eqn.  (6).  The  success  of  sediment  Fe  or  organic  carbon  as  the 
normalizing  factor  in  these  studies  could  be  due  to  its  dominance  in  the  sorption  of  the 
particular  trace  metal  investigated.  For  example,  the  strong  association  of  Pb  (Lion  et  al.  1982; 
Balistrieri  and  Murray  1982;  Leckie  et  al.  1984)  and  of  As  (Aggett  and  Roberts  1986;  Belzile  and 
Tessier  1990)  with  sedimentary  iron  oxyhydroxides  is  well  documented;  similarly,  the  affinity  of 
Hg  for  sediment  organic  matter  (NRCC  1988)  is  well  known. 
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From  the  preceding  discussion,  it  follows  that  the  observation  of  a  strong  relationship 
between  trace  metal  levels  in  a  benthic  species  and  sediment  characteristics  (e.g.,  extractable 
metal  concentration  normalized  with  respect  to  the  concentration  of  a  sediment  component  like 
Fe  oxyhydroxides  or  organic  matter)  cannot  be  considered  as  evidence  that  the  main  route  of 
trace  metal  uptake  is  via  ingestion  of  particulate  material.  Such  relationships  could  also  be 
observed  for  organisms  that  obtain  metal  directly  from  the  water,  provided  that  i)  the  pH  is 
approximately  constant  over  the  study  area;  ii)  the  dissolved  [M]  in  the  water  to  which  the 
organisms  are  exposed  is  in  sorptive  equilibrium  with  the  sediment  component  used  for 
normalization.  A  logical  consequence  is  also  that  the  ratios  {M}/{Fe}  or  {M}/{OM},  although 
much  better  predictors  than  total  sediment  metal  concentrations  in  some  cases,  are  not 
"universal"  predictors  in  the  sense  that  their  application  should  be  site-dependent  (i.e., 
restricted  to  a  narrow  pH  range). 


5.  CONCLUDING  REMARKS  AND  RECOMMENDATIONS 

We  strongly  believe  that  the  development  of  deterministic  models  bzised  on  sound  chemical 
and  biological  principles  constitutes  the  most  promising  approach  to  the  prediction  of  trace 
metal  burdens  in  benthic  organisms,  and  eventually  to  the  prediction  of  metal  effects  on  this 
community.  Provided  such  models  explicitly  incorporate  the  influence  of  major  environmental 
variables  on  metal  uptake,  they  should  be  of  general  applicability,  i.e.  they  should  be  able  to 
predict  the  metal  burden  of  benthos  in  lakes  other  than  those  used  for  their  calibration.  In 
contrast,  purely  empirical  models  (black  box  models),  based  on  statistical  relationships  only, 
are  useful  means  of  reducing  the  experimental  data  but  cannot  be  used  to  make  predictions 
beyond  the  original  data  set.  The  development  of  deterministic  models  has  been  hindered  in 
the  past  by  our  incomplete  knowledge  of  the  geochemical  and  biological  processes  that  control 
metal  accumulation  in  aquatic  organisms.  There  is  no  doubt  that  considerable  effort  should  be 
devoted  to  understanding  these  biogeochemical  processes  if  we  are  to  understand  how 
remedial  actions  will  affect  benthic  organisms. 

The  present  study  represents  an  attempts  to  develop  a  deterministic  model  of  general 
applicability  to  predict  bioaccumulation  of  trace  metals  in  aquatic  organisms.  In  the  formulation 
of  the  model,  we  use  concepts  derived  from  the  free-ion  model  of  metal-organism  interactions, 


and  from  surface  complexation  theory,  to  relate  cadmium  concentrations  in  bivalves  to  those  in 
the  water  or  in  the  surficial  oxic  sediments.  From  a  practical  point  of  view,  this  approach  has 
yielded  predictive  equations  for  Cd  accumulation  in  A.  grandis  that  require  as  input 
geochemical  variables  readily  measured  in  the  bottom  waters  (pH)  and  in  the  sediments 
({Cd}.j.'  {Fe-ox}  and  {OM}).  A  most  important  point  arising  from  the  study  is  that  relatively  little 
would  be  needed  to  establish  predictive  equations  for  metal  burdens  in  benthic  organisms  once 
the  geochemical  constants  x,  y,  Np^'*Kpj_^j  and  Nq^  • ''Kqj^,^^  are  known  accurately  for  the 
metals  of  interest. 

Several  recommendations  arise  from  our  work: 

1.  Determination  of  the  geochemical  constants.  The  approach  developed  here  for  Cd  could 
also  be  extended  to  other  trace  metals.  Since  the  geochemical  constants  x,  y,  Np^  •  *Kp^^j  and 
^OM  *  *^oM-cd  ^^  involved  in  the  prediction  of  cadmium  accumulation  by  (and  cadmium  effects 
on)  various  aquatic  organisms,  we  think  that  appreciable  effort  should  be  devoted  to  their 
accurate  determination;  such  an  effort  would  yield  a  basis  for  good  prediction  of  cadmium 
availability  to  benthic  organisms.  For  example,  the  following  points  are  suggested: 

-  Time-averaged  values  (rather  than  single  grab  measurements;  see  Appendix)  of  dissolved 
cadmium  concentrations  should  be  obtained  at  each  site  for  the  determination  of  the 
geochemical  constants. 

-  Possible  means  to  improve  the  estimation  of  the  free  aquo  ion  concentration,  [Cd^"^],  should 
be  investigated.  In  the  present  study,  we  have  considered  only  dissolved  inorganic  ligands;  the 
role  of  dissolved  organic  matter  in  complexing  Cd  should  be  assessed  and  taken  into  account . 

-  Our  knowledge  of  the  precise  sediment  components  responsible  for  the  sorption  of  Cd  should 
be  improved.  Which  forms  of  organic  matter  are  important  for  Cd  sorption?  In  the  present 
study,  we  assumed  that  it  is  total  organic  matter,  since  we  use  sedimentary  organic  carbon  to 
estimate  {OM}.  Which  forms  of  iron  oxyhydroxides  are  important  for  Cd  sorption?  The  reagent 
(NHjOH'HCL;  96  °C)  used  in  the  present  study  to  estimate  {Fe-ox}  extract  the  sum  of 
amorphous  iron  oxyhydroxides,  ferrihydrite  and  lepidocrocite.  More  precise  knowledge  of  the 
forms  of  organic  matter  and  iron  oxyhydroxides  responsible  for  the  sorption,  and  the 
development  of  appropriate  methods  to  measure  them,  would  yield  better  geochemical 
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constants.  Are  other  sediment  components  (e.g.  Mn  oxyhydroxides)  sometimes  important  for 
Cd  sorption?  In  the  present  study,  we  assumed  that  only  organic  matter  and  Fe  oxyhydroxides 
were  imp>ortant  for  sorbing  Cd. 

-  The  number  of  sites  where  Teflon  sheets  are  used  to  collect  diagenetic  iron  oxyhydroxides 
(and  associated  Cd)  should  be  increased;  these  experimental  data  are  needed  to  calculate  the 
geochemical  constants  x  and  N^^  «"Kp^^j,  and  hence  indirectly  to  calculate  y  and  Nq^  «'Kqj^^^ 
(see  Appendix  for  the  details  of  the  calculation).  In  the  present  study,  only  a  limited  number  of 
sites  (8)  were  investigated  with  the  Teflon  collection  technique. 

2.  Refinement  of  the  geochemical  model.  Note  that  the  considerations  raised  for  Cd  could 
also  be  extended  to  other  trace  metals.  The  prediction  of  [Cd(org)]  from  sedimentary 
measurements  depends  on  how  accurately  the  geochemical  model  can  relate  [Cd^"^]  in  the 
ambient  water  to  sedimentary  parameters.  The  geochemical  model  used  here  is  still  relatively 
crude,  and  further  work  to  refine  it  can  be  suggested. 

-  Our  understanding  of  the  nature  of  the  sorption  processes  could  be  improved.  Is  it  adsorption, 
surface  precipitation,  co-precipitation?  In  the  present  study,  we  assumed  that  [Cd^"*^]  is  related 
to  sediment  components  by  overall  surface  complexation  reactions  (see  Eqn  2  in  the  Appendix) 
suggested  for  adsorption  of  trace  metals  on  hydrous  metal  oxides  (Benjamin  and  Leckie,  1981). 
Identification  of  the  precise  processes  might  lead  to  better  expressions  relating  [Cd^"*"]  to 
sedimentary  components. 

-  Effort  should  be  devoted  to  identify  the  solid  phases  responsible  for  binding  [Cd].  Once  they 
are  known,  better  geochemical  models  could  be  developed  to  relate  [Cd^"^]  to  sediment 
components.  For  example,  all  the  experimental  data  on  the  adsorption  of  Cd  on  a  given  iron 
oxyhydroxide  can  be  described  with  a  set  of  intrinsic  acidity  constants,  intrinsic  adsorption 
constants  for  Cd  on  this  solid,  density  of  sites  and  specific  surface  area  of  this  solid,  and  a 
description  of  the  potential-surface  charge  relationship  (Dzombak  and  Morel,  1990).  Laboratory 
experiments  with  pure  solid  phases  have  shown  that  values  of  these  parameters  vary  from  one 
Fe  oxyhydroxide  to  the  other.  If  the  exact  forms  of  Fe  oxyhydroxides  present  in  natural 
sediments  were  known,  we  could  use  these  laboratory-derived  parameters  in  a  refined 
geochemical  model.  Since  this  information  is  not  yet  available,  we  adopted  in  the  present  study 
a  semi-empirical  approach  where  an  overall  reaction  (see  Eqn  2  in  the  Appendix)  is  used  to 
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describe  sorption  of  Cd  on  natural  Fe  oxy hydroxides. 

3.  Biomonitors.  Aquatic  organisms  have  been  proposed  as  "biomonitors"  for  trace  metal 
pollution.  Our  results  raise  several  points  and  recommendations  for  further  work. 

-  Freshwater  bivalves  have  been  proposed  for  biomonitoring  purposes  (Forester,  1980). 
According  to  our  results,  A.  grandis,  and  probably  E.  complanata,  do  not  react  rapidly  to 
changes  in  ambient  Cd,  Cu  and  Pb  concentrations.  For  example,  A.  grandis  specimens 
transferred  to  a  polluted  lake  had  not  attained,  after  three  years,  the  metal  burdens  of 
indigenous  specimens  in  the  polluted  lake.  Thus,  freshwater  bivalves  will  probably  be  useful  for 
monitoring  long-term  trends  in  metal  pollution.  The  suitability  of  smaller  organisms  (e.g.  insect 
larvae)  should  be  investigated  for  monitoring  trace  metal  variations  in  the  environment  on  a 
short-term  basis. 

-  Use  of  organisms  for  biomonitoring  trace  metals  requires  that  the  relationships  between  metal 
concentrations  in  the  organisms  (or  in  a  specific  tissue)  and  their  external  environment:  i)  be 
well  known;  ii)  apply  to  sites  other  than  those  used  for  calibration;  iii)  apply  to  environmental 
conditions  other  than  those  used  for  the  calibration.  In  our  opinion,  this  generality  can  be 
obtained  only  if  the  relationships  between  tissue  metal  concentration  and  predictive 
environmental  variables  are  derived  from  sound  biological  and  chemical  principles.  This  critical 
point  has  been  overlooked  in  the  past  and  should  be  addressed  explicitly. 

-  Figure  3  shows  that  prediction  of  Cd  concentration  in  the  tissues  oiA.  grandis  from 
sedimentary  variables  is  acceptable  for  a  large  number  of  aquatic  environments  (over  a 
geographical  area  of  about  350  0(X)  km^  in  this  case);  this  is  an  example  of  a  predictive 
relationship  based  on  biological  and  chemical  processes.  Admittedly,  the  relationships  between 
[Cd(org)]  and  the  sedimentary  parameters  could  be  improved  by  a  better  determination  of  the 
geochemical  constants  and  a  refinement  of  the  geochemical  model,  as  discussed  above. 
Similar  approaches  to  the  one  developed  for  the  study  of  Cd  in  /I.  grandis  could  be  adopted  for 
the  study  of  other  trace  metals  in  other  aquatic  organisms. 
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Predicting  metal  bioaccumulation  in  ben- 
thic  organisms  as  a  function  of  environmental 
variables  is  not  straightforward.  Many  of  these 
organisms  are  in  contact  with  both  dissolved 
and  particulate  trace  metals  and  can  in  prin- 
ciple accumulate  the  metals  either  directly  from 
the  water  or  from  food  particles  (Luoma  1 983). 
Uptake  from  either  source  will  be  influenced 
by  physicochemical  factors  in  the  aqueous  (e.g. 
pH,  ligand  type,  and  concentrations)  and  par- 
ticulate phases  (e.g.  association  of  the  trace 
metal  with  specific  sediment  constituents).  The 
issue  is  further  complicated  by  the  fact  that 
dissolved  metal  concentrations  in  the  bottom 
waters  are  related  to  those  present  in  the  un- 
derlying surficial  sediments  (Tessier  1992). 

Considerable  uncertainties  remain  in  our 
knowledge  of  the  processes  responsible  for  the 
regulation  and  cycling  of  trace  metals  in  lakes. 
However,  comparison  of  water-column  pro- 
files of  trace  metals  and  nutrients  (Murray 
1987)  and  chemical  analysis  of  material  col- 
lected with  sediment  traps  (Sigg  1987)  suggest 
that  trace  metals  (Cd,  Cu,  Pb,  Zn)  can  be  re- 
moved efficiently  by  phytoplankton  in  circum- 
neutral  lakes  and  transported  to  the  bottom 
sediments  by  the  sinking  organisms.  Sediment 
trap  measurements  suggest  that  Fe  and  Mn 
oxyhydroxides  are  also  important  scavengers 
of  trace  metals  (Cd,  Cu,  Pb,  Zn;  Sigg  1987; 
Sigg  et  al.  1987)  and  can  carry  them  to  the 
bottom  sediments  in  lakes.  Transport  of  trace 
metals  to  lake  sediments  is  not  limited  to  par- 
ticulate forms;  dissolved  metals  can  also  dif- 
fuse from  the  bottom  waters  into  the  sediments 
in  acid  lakes  (Cu,  Ni,  Zn;  Carignan  and  Nriagu 
1985;  Carignan  and  Tessier  1985;  Tessier  et 
al.  1989).  Release  of  metals  from  the  sediments 
back  to  the  overlying  waters,  following  redox 
processes,  may  also  occur  (Morfett  et  al.  1988). 
Thus,  trace  metals  in  surficial  oxic  sediments 
are  expected  to  be  associated  with  sediment 
components  such  as  organic  matter  (living  or- 
ganisms, debris,  humic  acids)  and  Fe  and  Mn 
oxyhydroxides,  and  the  trace  metal  content  of 
these  sediment  components  can  be  expected 
to  be  related  to  that  in  the  overlying  water. 

The  above  considerations  stress  the  impor- 
tance of  understanding  geochemical  processes 
if  general  models  for  predicting  metal  accu- 
mulation in  benthic  organisms  are  to  be  de- 
veloped. In  this  paper  we  illustrate  the  point, 
taking  as  an  example  the  accumulation  of  Cd, 


a  nonessential  trace  metal,  in  the  freshwater 
pelecypod  Anodonta  grandis,  a  benthic  filter 
feeder,  for  a  large  number  of  lakes  over  a  wide 
geographical  area.  We  discuss  the  relationships 
between  concentrations  of  Cd  in  oxic  sedi- 
ments and  those  in  the  overlying  water,  by 
means  of  surface  complexation  concepts.  We 
then  use  these  concepts  together  with  the  free- 
metal  ion  model  (Morel  1983)  to  relate  Cd 
concentrations  in  A.  grandis  to  dissolved  and 
sedimentary  Cd  concentrations  in  the  animals' 
environment. 

Study  area 

Oxic  sediments,  associated  pore  water,  and 
bivalves  (when  present)  were  collected  at  49 
littoral  stations  in  38  lakes  distributed  over  a 
350,000-km-  study  area  in  Québec  and  On- 
tario (Table  1).  The  sites  were  chosen  to  rep- 
resent a  gradient  of  Cd  contamination  (Cd  con- 
centrations in  water,  sediments,  and  bivalves) 
and  a  wide  range  of  lake  pH  values  and  other 
chemical  characteristics  (Table  2). 

For  the  lakes  in  the  Eastern  Townships  of 
Québec,  the  bedrock  is  composed  of  Ordovi- 
cian  and  Cambrian  sedimentary  rock,  covered 
by  unconsolidated  deposits  of  glacial  till  and 
marine  clays;  small-scale  mining  activities  (Cu, 
Ni,  Pb,  Zn)  occurred  until  recently.  All  the 
other  lakes  are  located  on  the  Precambrian 
Shield.  Lakes  in  the  mining  areas  of  Sudbury 
(Ni,  Cu)  and  Rouyn-Noranda  (Cu,  Pb,  Zn)  are 
subject  to  relatively  high  acid  and  metal  con- 
tamination from  current  mining  operations, 
abandoned  mines,  and  atmospheric  deposi- 
tion from  nearby  smelters.  Lakes  in  the  area 
of  Chibougamau  are  also  affected  by  ongoing 
mining  operations  (Cu).  Lakes  in  the  Muskoka 
area  are  far  from  industrial  activities  (~200 
km  southeast  of  Sudbury  and  250  km  north 
of  Toronto);  Lake  Tanlaré  is  40  km  north  of 
Québec  city,  entirely  within  an  ecological  re- 
serve. 

Materials  and  methods 

Sampling— 'Waier  samples  were  collected  in 
the  first  few  centimeters  above  the  sediment- 
water  interface  with  in  situ  samplers  (pore- 
water  peepers;  1-cm  vertical  resolution;  Gel- 
man  HT-200  membrane:  three  per  site)  similar 
to  those  described  by  Hesslein  (1976)  and  Ca- 
rignan et  al.  (1985).  The  peepers  consist  of 
Plexiglas  sheets  (1.3  cm  thick)  into  which  com- 
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Sample  identification  and  location  of  stations. 


A-OI 

Aylmer 

Eastern  Townships,  Qu 

A-02 

BC-OI 

Blue  Chalk 

Muskoka,  Onlano 

BE-03 

Beauchastel 

Rouyn-Noranda,  Québ< 

Bl-OI 

Bird 

Muskoka 

BO-o; 

Bousquet 

Rouyn-Noranda 

BR-01 

La  Bruère 

Rouyn-Noranda 

BR-0-1 

BR-05 

BW-01 

Bigwind 

Muskoka 

BW-02 

CA-01 

Caron 

Rouyn-Noranda 

CE-05 

Brompton 

Eastern  Townships 

CE-06 

' 

CH-01 

Chub 

Muskoka 

CL-03 

Clearwater 

Sudbury.  Ontario 

DA-01 

D'Alemben 

Rouyn-Noranda 

DF-01 

Dufresnoy 

Rouyn-Noranda 

DU-01 

Dufay 

Rouyn-Noranda 

FA-02 

Fairbank 

Sudbury 

FL-01 

Flavnan 

Rouyn-Noranda 

GF-01 

Gullfeaiher 

Muskoka 

GF-02 

H-01 

Harp 

Muskoka 

HE-OI 

Héva 

Rouyn-Noranda 

J-01 

Joannes 

Rouyn-Noranda 

J-02 

LB-01 

Lake  of  Bays 

Muskoka 

LD-OI 

Lac  aux  Dorés 

Chibougamau,  Québec 

LD-02 

LD-03 

LD-04 

MA-01 

Massawippt 

Eastern  Townships 

MC-02 

McFarlane 

Sudbury 

ME-OI 

Memphrémagog 

Eastern  Townships 

NO-01 

Norhan 

Chibougamau 

OP-01 

Opasatica 

Rouyn-Noranda 

PA-01 

Paint 

Muskoka 

RI-01 

Ril 

Muskoka 

SC-01 

St.  Charles 

Sudbury 

SI-01 

Silver 

Sudbury 

SN-OI 

St.  Nora 

Muskoka 

TA-03 

Taniaré 

Québec,  Québec 

TA-04 

TI-01 

Tilton 

Sudbury 

TO-01 

Tock 

Muskoka 

VA-01 

Vaudray 

Rouyn-Noranda 

WA-01 

Wavy 

Sudbury 

WR-01 

Little  Wren 

Muskoka 

45°50'N,  71''20'W 

45''I2'N,  78°56W 
48*I0'N,  79t)8'W 
45°02'N,  79'D4'W 
48''13'N,  78''39'W 
48°09'N,  78*56'W 


45°03'N,  78°50'W 


47°55'N,  78°58'W 
45°25'N,  72°I0'W 


45''13'N, 
46''22'N. 
48r23'N, 
48"'27'N, 
48't)3'N, 
46'>27'N. 
48*18'N, 
45°06'N, 


78''59'W 
8I"03'W 
79n)l'W 
79°02'W 
79°28'W 
8r25'W 
79''12'W 
79't)l'W 


45'23'N,  79't)8'W 
48*'irN,  78''19'W 
48"'I2'N,  78''40'W 


45''15'N, 
49'55'N, 
49*49'N, 
49°55'N, 
49*54'N, 
45°15'N. 
46°25'N, 
45°05'N, 
49°29'N, 
48°05'N, 
45''13'N, 
45°10'N, 
46°26'N, 
46°25'N, 
45''10'N, 
47»04'N, 


79°00'W 
74°I6'W 
74'20'W 
74°14'W 
74''I7'W 
72'^0'W 
80°57'W 
72''IO'W 
74°26'W 
79°18'W 
78''57'W 
79°O0'W 
8I°01'W 

8i°o.rw 

78'50'W 
71»32'W 


46'*22'N,  8I°04'W 
45°16'N,  78''53'W 
48°07'N,  78''42'W 
46°17N,  81°06'W 
45''11'N,  78-5  rw 


partments  (3.3  ml,  1  cm  apart  in  a  row,  or- 
ganized in  two  parallel  rows)  have  been  ma- 
chined. When  in  storage,  the  peepers  were  kept 
in  a  dilute  H^SOj  solution.  Before  use,  they 
were  nnsed  with  demineralized  water  (Milli- 
Q3RO/Milli-Q2  system;  Millipore  Ltd.);  the 
compartments  were  filled  with  demineralized 


water  and  deoxygenaled  by  bubbling  with  ni- 
trogen for  at  least  24  h  in  Plexiglas  cylinders 
filled  with  demineralized  water  before  being 
inserted  vertically  in  the  lake  sediments  by 
SCUBA  divers.  After  a  2-week  equilibration 
period  in  the  sediments,  the  peepers  were  re- 
trieved by  divers  and  sampled  immediately. 
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Table  2.  Mean  concentrations  of  total  dissolved  Cd  ([Cd)).  free  Cd  ([Cd-*]),  major  ions,  alkalinity,  and  organic  C 
((C„^l)  in  the  overlying  waters,  and  concentrations  of  total  Cd  (iCdlj),  Cd  associated  with  organic  matter  ((OM-Cdl; 
calculated  with  Eq.  16),  Fe  oxyhydroxides  ((Fe-ox)),  and  organic  matter  ((OM);  Mmol  Cg"')  in  the  surficial  sediments 
at  the  various  sites  studied.  (Not  measured  — NM.) 


ICa) 


|Na| 


[SO.l         (Alkl* 


IC„1 


ICdl 


lCd-1 


iCdl, 


lOM-Cdl  lFff-o«l 


pH 


A-OI 

A-02 

BC-OI 

BE-03 

BI-01 

BO-02 

BR-01 

BR-04 

BR-05 

BW-01 

BW-02 

CA-OI 

CE-05 

CE-06 

CH-01 

CL-03 

DA-OI 

DF-01 

DU-01 

FA-02 

FL-OI 

GF-01 

GF-02 

H-OI 

HE-01 

J-01 

J-02 

LB-OI 

MA-OI 

MC-02 

ME-OI 

NO-01 

OP-01 

PA-01 

RI-01 

SC-01 

SI-01 

SN-01 

TA-03 

TA-04 

TI-01 

TO-OI 

VA-OI 

WA-01 

WR-OI 

LD-01 

LD-02 

LD-03 

LD-04 


7.25 
7.11 
6.82 
7.63 
6.37 
6.34 
7.94 
7.33 
7.33 
6.05 
6.52 
7.05 
7.16 
7.16 
5.50 
4.79 
7.14 
7.17 
6.55 
7.36 
7.39 
5.90 
6.07 
6.93 
6.19 
7.27 
7.27 
6.94 
8.10 
7.49 
7.56 
7.09 
7.39 
7.13 
6.02 
6.88 
4.07 
6.80 
5.68 
5.65 
5.88 
6.43 
6.57 
4.62 
6.27 
7.42 
7.50 
7.47 
7.40 


124 

128 

57 

484 

107 

106 

395 

393 

369 

94 

86 

245 

102 

100 

47 

136 

197 

166 

85 

163 

158 

84 

67 

59 

52 

176 

188 

64 

700 

410 

380 

106 

280 

147 

113 

211 

218 

57 

NM 

NM 

97 

139 

90 

52 

80 

176 

224 

177 

200 


74 
72 
26 

NM 
35 
45 

NM 

NM 

248 
46 
52 
105 

NM 
39 
22 
94 
62 

NM 

NM 
38 

NM 
43 
20 
41 
40 

NM 
35 
44 
180 
1,435 
165 
25 
111 
86 
44 

758 

1,780 

54 

NM 

NM 
55 
39 

NM 
28 
64 
30 
48 
29 
39 


57 

64 

65 

NM 

78 

82 

NM 

NM 

241 

72 

75 

181 

NM 

62 

71 

175 

108 

NM 

NM 

144 

NM 

81 

86 

67 

68 

NM 

101 

77 

108 

227 

81 

31 

120 

97 

78 

197 

439 

66 

48 

50 

152 

79 

NM 

126 

81 

47 

92 

48 

69 


507 
453 
135 
680 
109 

89 
579 
687 
570 
145 

90 
276 
364 
581 

32 
0 
267 
NM 
NM 
246 
NM 

79 

41 
175 

40 
NM 
679 
128 
1,430 
669 
882 
178 
494 
112 

57 
187 
0 
111 
NM 
NM 

37 

61 

NM 

0 

50 
377 
401 
418 
404 


767 
758 
363 
NM 
475 

1,217 
NM 
NM 
950 
524 
625 

1,108 
NM 
453 
618 
343 

1,003 
NM 
NM 
375 
NM 
858 
787 
593 

1,000 
NM 
784 
470 
398 
482 
401 
813 
818 
440 
473 
645 
425 
463 
NM 
NM 
501 
474 
NM 
420 
453 
462 
492 
448 
524 


0.110 
0.089 
0.041 
NM 
0.170 
2.37 
NM 
NM 
0.790 
0.350 
0.21 
3.14 
NM 
0.098 
0.320 
4.83 
1.77 
NM 
NM 
0.036 
NM 
0.21 
0.10 
0.140 
0.67 
NM 
1.13 
0.090 
0.049 
0.336 
0.032 
0.114 
0.37 
0.043 
0.20 
0.37 
17.5 
0.071 
0.537 
0.400 
1.93 
0.23 
NM 
2.32 
0.25 
0.018 
0.062 
0.027 
0.024 


0.108 
0.087 
0.041 

0.167 

2.32 


0.743 
0.344 
0.206 
3.01 

0.096 
0.314 
4.55 
1.72 


0.035 

0.21 
1.10 
0.138 
0.66 

0.10 
0.088 
0.047 
0.285 
0.031 
0.113 
0.36 
0.042 
0.20 
0.33 
14.1 
0.069 
0.537 
0.400 
1.85 
0.22 

2.16 

0.244 

0.018 

0.060 

0.027 

0.023 


4.00 
3.60 
2.80 
49.8 

8.43 
33.0 
52.1 
43.6 
36.1 
5.00 
2.33 
51.2 
10.5 
8.29 
1.90 
3.87 
93.3 
18.3 
9.13 
3.27 
27.8 
4.19 
8.91 
9.00 
10.4 
62.8 
61.9 
6.00 
3.40 
108 
6.60 
9.92 
4.40 
2.67 
9.79 
52.5 
1.34 
5.40 
5.27 
7.15 
8.78 
18.0 
38.3 
3.31 
11.2 
6.92 
9.32 
3.10 
15.0 


3.64 
3.41 

2.52 

S.09 
30.1 


33.6 
4.86 
2.26 

48.0 

8.09 
1.67 
2.89 
91.3 
17.5 
8.24 
3.07 
27.0 
3.98 
8.58 
8.72 
9.35 

60.1 
5.59 
3.34 
103 
6.43 
9.82 
3.76 
2.38 
9.37 

50.2 
1.10 
5.00 
3.42 
6.64 
6.93 
17.6 
34.9 
2.68 

II. I 
6.59 
9.13 
3.04 

14.9 


157 
26.9 
65.6 

143 

102 

207 
94.2 
75.1 
72.9 
30.6 
24.8 
92.9 

103 
97.7 
57.7 

757 

119 

126 
88.0 
61.1 

123 
46.7 
25.9 
27.5 

122 
94.9 
.  120 
42.6 

117 

371 

167 
98.2 
40.8 

112 
83.7 

233 

544 

115 
3,212 

482 

590 

120 

129 

683 
71,7 

165 

139 
43.7 
88.5 


1,775 

575 
■  775 

NM 
4,472 
3,413 

NM 

NM 
1,093 
1,804 
1,200 
1.707 

NM 
4,617 

883 
5,967 
6,533 
3.251 
1,175 
1,008 
4,198 
1.601 
1,159 
1.075 
1,838 

NM 
4,675 

733 
5.592 
7,883 
6.642 
12,158 

262 
1.107 
3,230 
6,575 
8,375 
1,933 
11.633 
1 2.400 
4,042 
9,967 
1.883 
8.275 
9.358 
3.550 
6,883 
2,458 
16.125 


•  licq  In 


ulalcd  wiih  HYDRAQL,  given  Ihc  inorganic  C 


and  ihc  other  chemical  charactenstics  of  the  lakei 


Samples  (1  ml)  in  the  first  row  were  collected 
from  five  compartments  above  the  water-sed- 
iment interface,  and  pH  was  measured  in  each 
sample  in  the  field  (Carignan    1984)  with  a 


combined  microelectrode  (Microelectrodes, 
Inc.,  model  MI-710)  and  a  portable  pH  meter 
(Radiometer,  model  PHM80).  On  several  oc- 
casionsdunng  the  summer  season,  water  sam- 
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pies  were  also  collected  at  many  of  the  sites 
(33/49)  with  a  clean  polyethylene  bottle  close 
to  the  sediment-water  interlace  for  pH  mea- 
surement. The  purpose  of  this  measurement 
was  to  calculate  time-averaged  pH  values 
[  — logio(2[H"l/Ai))  for  the  summer  period. 
Samples  (I  ml)  for  dissolved  sulfate  and  chlo- 
ride analysis  were  removed  from  the  peepers 
with  a  syringe  from  compartments  of  the  same 
row  and  injected  into  prewashed  polypropyl- 
ene tubes;  those  for  inorganic  C  determina- 
tions (1  ml)  were  also  obtained  from  the  same 
row  with  a  syringe  and  injected  through  a  sep- 
tum into  pre-evacuaied  and  prewashed  glass 
tubes. 

The  samples  (3.3  ml)  for  trace  metal  analysis 
were  collected  from  the  compartments  of  the 
second  row  by  piercing  the  peeper  membrane 
with  a  Gilson  pipette  fitted  wdth  an  acid-cleaned 
tip:  these  samples  were  injected  into  pre- 
washed and  preacidified  (30  m1  of  1  N  Ulirex 
HNOj,  final  pH  < 2.5)  Teflon  vials.  Additional 
water  was  collected  from  the  compartments  6 
and  7  cm  above  the  sediment-water  interface 
of  each  peeper  (combined  content,  5  ml)  for 
analysis,  if  necessary,  by  multiple  injection 
flameless  atomic  absorption  spectrophotom- 
eti7  (AAS,  see  below).  Mean  concentrations  of 
total  dissolved  Cd  ([Cd])  in  the  overlying  wa- 
ters are  given  for  each  site  in  Table  2.  When 
sufficiently  high  (  >  0. 7  nM),  the  concentrations 
found  in  the  five  compartments  of  the  pore- 
water  peepers  immediately  above  the  sedi- 
ment-water interface  (i.e.  3x5=15  values) 
were  used  for  calculating  the  means;  in  the 
other  cases,  the  values  obtained  by  successive 
injections  (i.e.  3  values)  were  used  for  the  cal- 
culation. Inadvertent  contamination  of  the 
compartments  with  Cd,  occurring  before  in- 
serting the  peepers  into  the  sediments  despite 
the  precautions  taken,  would  have  decreased 
by  diffusion  out  of  the  cells  during  the  2-week 
equilibration  period  in  the  lake  sediments.  A 
greater  possibility  of  contamination  with  these 
devices  probably  occurs  during  sampling  of  the 
compartments,  since  small  particles  (contain- 
ing Cd)  that  adhere  to  the  membrane  can  be 
taken  inadvertently  along  with  the  water  sam- 
ple. The  low  values  of  [Cd]  found  m  many 
circumneutral  pH  lakes  (10""  M  range;  see 
Table  2)  argue  against  important  contamina- 
tion problems. 

Sediment  cores  were  collected  bv  divers. 


close  to  the  peepers,  with  Plexiglas  tubes  (9- 
cm  diam).  The  tubes  were  tightly  closed  to 
minimize  perturbation  of  the  sediments  dunng 
their  transport  to  the  shore.  The  sediment  cores 
were  extruded  on  shore,  and  only  the  upper- 
most 0.5  cm,  containing  oxidized  sediments, 
was  retained.  These  samples  were  placed  in 
500-ml  centrifugation  bottles  half  filled  with 
lake  water  and  kept  at  ~4°C  during  transport 
to  the  laboratory  where  they  were  kept  frozen 
until  analysis. 

In  several  of  the  lakes  (sites  CE-05,  CL-03, 
DA-0 1 ,  HE-0 1 ,  J-02,  MC-02,  SI-0 1 ,  and  WA- 
01),  Teflon  sheets  were  inserted  venically  in 
the  sediments  by  divers  and  left  in  place  for 
more  than  10  weeks.  In  sufficiently  reducing 
sediment  layers,  Fe(III)  oxyhydroxides  are  re- 
duced to  Fe(II)  (due  to  their  burial  or  to  fluc- 
tuation of  the  redox  transition  zone),  which  is 
released  to  the  fwre  water;  a  portion  of  the 
Fe(II)  thus  produced  diffuses  upward  where  it 
is  oxidized  in  the  upper  sediment  layers  and 
hydrolyzes  to  form  diagenetic  Fe(III)  oxyhy- 
droxides. The  diagenetic  Fe(III)  oxyhydrox- 
ides thus  produced  can  be  collected,  along  with 
the  sorbed  trace  elements,  on  inert  material 
such  as  sheets  of  Teflon  (Belzile  et  al.  1989; 
De  Vitre  et  al.  1991). 

Fe  oxyhydroxide  samples  obtained  by  this 
in  situ  technique  appear  as  a  thin  orange-brown 
band,  typically  0.5  cm  wide  and  5-20  nm  thick 
for  littoral  sites  (Fortin  et  al.  in  prep.).  The 
amount  of  material  deposited  is  small  (typi- 
cally 500  Mg  per  Teflon  sheet),  but  it  has  the 
advantage  of  being  relatively  free  from  con- 
tamination by  the  sediment  matrix.  The  Fe 
particles  collected  have  been  characterized  by 
transmission  electron  microscopy  and  by  elec- 
tron and  X-ray  diffraction.  Their  chemical 
composition  has  been  determined  at  a  micro- 
scopic scale  by  energy-dispersive  spectroscopy 
and  on  a  macroscopic  scale  by  wet  chemical 
methods  (Fortin  et  al.  in  prep.).  According  to 
these  analyses,  the  crystallinity  of  the  panicles 
is  low;  crystalline  forms  identified  are  poorly 
ordered  ferrihydrite  and  lepidocrocite.  Si,  SO4, 
CI,  phosphate,  Mn,  Ca,  and  Al  represent  minor 
components  of  the  Fe-rich  particles,  whereas 
organic  C  (presumably  microorganisms,  as 
suggested  from  transmission  electron  micros- 
copy observations)  is  a  more  abundant  com- 
ponent. After  retrieval  of  the  Teflon  collectors, 
the  areas  covered  with  the  Fe  oxyhydroxides 
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were  cut  out  from  the  Teflon  sheet  and  the 
deposited  material  was  dissolved  at  96°C  in  a 
0.04  M  NHjOHHCl  solution  in  25%  (vol/ 
vol)  acetic  acid.  Fe  and  Cd  concentrations  in 
these  extracts  were  determined  by  AAS  as  de- 
scribed below  for  the  sediment  extracts. 

Specimens  of  A.  grandis  (usually  10  per  site; 
generally  ranging  from  8  to  10  cm  long)  were 
obtained  by  divers  at  each  site,  within  a  radius 
of  50  m  from  the  sediment  and  water  collection 
site.  The  bivalves  were  placed  in  plastic  bags 
with  lake  water  and  maintained  at  ~4°C  dur- 
ing transport  to  the  laboratory  where  they  were 
left  to  depurate  for  at  least  24  h  in  aerated  lake 
water.  They  were  then  dissected  into  gills, 
mantle,  hepatopancreas,  and  remaining  tis- 
sues, hereafter  referred  to  as  "remains."  No 
organisms  were  gravid,  as  sampling  was  un- 
dertaken early  in  summer.  For  each  sampling 
site,  each  tissue  type  of  the  10  animals  col- 
lected was  pooled  and  frozen  at  —  20°C  until 
needed  for  analysis. 

Analyses— The  water  samples  were  analyzed 
for  Cd,  Fe,  Mn,  SO4,  CI,  major  cations,  and 
organic  and  inorganic  C.  The  metal  concen- 
trations were  obtained  by  flame  AAS  when 
possible  (Fe,  Mn,  Ca,  Mg,  Na,  K;  Varian  Tech- 
tron  model  575ABQ  or  model  Spectra  AA-20) 
or,  otherwise,  by  flameless  AAS  (Fe,  Mn,  Cd; 
Varian  Techtron  model  1275  or  Spectra  AA- 
30;  GTA-95  or  GTA-96).  When  initially  mea- 
sured levels  of  Cd  in  the  overlying  waters  did 
not  exceed  a  concentration  threshold  of  five 
times  the  analytical  detection  limit  (0. 1 3  nmol 
liter"'),  its  concentration  was  determined  by 
a  multiple  injection  technique.  The  subsam- 
ples  were  injected  successively  and  subjected 
to  the  preliminary  drying  steps;  the  combined 
sample  was  then  atomized.  SO4  and  CI  con- 
centrations were  determined  by  ion  chroma- 
tography (Dionex  Autolon,  system  12);  dis- 
solved inorganic  C  was  measured  by  gas 
chromatography  (Carignan  1984)  and  dis- 
solved organic  C  (DOC)  by  persulfate-UV  ox- 
idation, followed  by  conductometric  deter- 
mination on  a  Technicon  AutoAnalyzer  of  the 
CO2  released.  The  equilibnum  model  HY- 
DRAQL  (Papelis  et  al.  1988)  was  used  to  cal- 
culate ion  activities.  The  calculation  of  Cd^* 
activity  involved  the  mean  [Cd]  in  the  over- 
lying water  (Table  2),  the  measured  concen- 
trations of  the  inorganic  ligands  OH",  COj^", 
SOj'^  and  CI   ,  and  the  stability  constants  of 


the  inorganic  complexes  (Smith  and  Martell 
1977). 

The  surficial  sediment  samples  were  thawed 
and  centrifuged  to  remove  excess  water;  sub- 
samples  (equivalent  to  ~  1  g  dry  wt)  were  ex- 
tracted, and  Cd  (and  Fe)  was  partitioned  into 
the  following  empirical  fractions  (Tessier  et  al. 
1989):  (1)  the  sediment  subsample  was  ex- 
tracted with  MgClj;  (2)  the  residue  from  (I) 
was  extracted  with  an  acetate  buffer  at  pH  5; 
(3)  the  residue  from  (2)  was  extracted  at  room 
temperature  with  NHjOHHCl;  (4)  the  resi- 
due from  (3)  was  extracted  with  NH^OHHCl 
at  96°C;  (5)  the  residue  from  (4)  was  extracted 
with  H2O2;  (6)  the  residue  from  (5)  was  di- 
gested with  a  mixture  of  hydrofluoric,  nitric, 
and  perchloric  acids.  Fe,  Mn,  and  Cd  concen- 
trations in  the  extracts  were  determined  by 
flame  AAS  with  the  appropriate  extractant 
matrixes  for  standards  and  blanks.  Details  of 
these  procedures  are  given  elsewhere  (Tessier 
et  al.  1979,  1989).  Sediment  organic  C  con- 
centrations, {Co,^},  were  determined  with  a 
CNS  analyzer  (Cario-Erba,  model  NA1500) 
after  removal  of  inorganic  C  by  acidification 
with  H2SO4  (0.5  mol  liter"',  15  min,  100  ml 
g"'  sediment  dry  wt). 

Each  pooled  bivalve  tissue  was  homoge- 
nized (Brinkman  tissue  grinder,  model  CH- 
60 1 0).  A  subsample  was  then  dried  to  constant 
weight  to  determine  the  wet  :  dry  wt  ratio,  and 
a  second  subsample  was  digested  in  a  Teflon 
bomb  with  concentrated  nitric  acid  (Aristar;  3 
ml  per  1 00  mg  of  tissue  dry  wt)  in  a  microwave 
oven  at  pressures  between  5,500  and  7,000  kPa 
for  ~  1  min.  Cd  concentrations  in  the  diluted 
digested  samples  were  determined  by  flame 
AAS.  A  certified  reference  material  (lobster  he- 
patopancreas, TORT-1,  Natl.  Res.  Council 
Canada)  was  regularly  submitted  to  the  same 
digestion  procedure  and  analyzed  for  Cd.  We 
measured  25.7±1.9  ^g  Cd  g"'  («  =  10)  for 
TORT-1  (certified  value,  26.3  ±2.1  Mg  g"'). 
Unless  otherwise  indicated,  reported  variabil- 
ities refer  to  standard  error. 

Results  and  discussion 

Relationship  between  dissolved  and  surficial 
sedimentary  Cd— Table  2  shows  that  the  mean 
[Cd]  varies  by  nearly  three  orders  of  magnitude 
over  the  study  area.  Lakes  in  the  areas  of  Rou- 
yn-Noranda  and  Sudbury  show  the  highest 
[Cd],  a  probable  consequence  of  the  smelting 
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operations  in  these  areas.  There  is,  however. 
a  general  tendency  for  [Cd]  to  decrease  with 
increasing  lake  pH  (Fig.  1  A),  in  the  absence  of 
local  Cd  sources.  For  example,  the  linear  re-         _ 
gression  between  log[Cd]  and  pH  is  '^ 

loglCdl=  -0.51(±0.07)pH  -  6.6(±0.2)(1)         ^ 


with  r  =  0.1\{n=  26)  if  lakes  from  the  Rouyn- 
Noranda  and  Sudbury  areas,  for  which  the 
smelters  constitute  a  known  anthropogenic 
source  of  Cd,  are  removed. 

For  all  of  the  lakes  studied  (pH  4.0-8.4;  al- 
kalinity concn  between  0  and  1,400  meq  li- 
ter"'; [PO4]  <  0. 1  mM;  [Si]  <  20  mM),  CdCOj(s) 
should  be  the  most  stable  Cd  solid  phase.  Cal- 
culation of  the  saturation  index  [SI  =  log(IAP/ 
A:J,  where  lAP  is  the  ion  activity  product  and 
K,  the  solubility  product]  indicates  a  clear  un- 
dersaturation  of  the  overlying  waters  of  all  the 
lakes  with  respiect  to  this  solid  phase  (Fig.  IB). 
In  general,  the  [Cd-^  ]  calculated  by  taking  into 
account  the  inorganic  ligands  represents  >95% 
of  [Cd]  (Table  2). 

Total  Cd  concentrations  in  the  sampled  sed- 
iments vary  by  more  than  two  orders  of  mag- 
nitude (ICd},;  Table  2);  the  highest  concen- 
trations are  found  in  lakes  in  the  mining  areas 
of  Rouyn-Noranda  and  Sudbury.  Because  sol- 
ubility equilibrium  involving  pure  solid  phases 
cannot  explain  the  presence  of  this  metal  in 
oxic  sediments,  other  reactions  including  ad- 
sorption, absorption,  surface  precipitation,  and 
coprecipitation  must  be  invoked  to  relate  its 
concentrations  in  the  overlying  water  to  its 
concentrations  in  the  surficial  sediment  (Tes- 
sier  1992).  The  general  term  "sorption"  is  used 
herein,  because  the  different  reactions  cannot 
be  distinguished  in  natural  waters  (Honeyman 
andSantschi  1988). 

There  are  several  indications  in  the  litera- 
ture that  Fe  and  Mn  oxyhydroxides  and  or- 
ganic matter  are  important  components  for 
sorbins  Cd  in  sediments.  For  example,  Luoma 
and  Bryan  (1981)  have  shown,  by  statistical 
analysis  of  50  oxidized  surficial  estuarine  sed- 
iments presenting  a  wide  range  of  sedimentary 
Cd  and  sediment  component  concentrations, 
that  sedimentary  Cd  was  correlated  with  ex- 
tractable  Fe  (0.2  M  ammonium  oxalate  in  0.2 
M  oxalic  acid  or  1  M  HCl);  extractable  Fe 
(presumably  Fe  oxyhydroxides)  was  found  to 
be  more  important  than  total  Fe  for  binding 
Cd.  Lion  et  al.  (1982)  examined  the  sorption 
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Fig.  1 .  (Cd)  and  saturation  index  as  a  function  of  pH 
for  the  various  sites  studied.  Different  symbols  are  used 
for  the  lakes  located  in  the  areas  of  Sudbury  (O),  Rouyn- 
Noranda  (•),  Chibougamau  (D),  Muskoka  (■),  Eastern 
Townships  (A),  and  Québec  (A).  Saturation  index  is 
log(IAP/A:,).  Dashed  line  indicates  saturation  with  respect 
to  CdCO,(s). 


of  Cd  on  oxidized  surficial  sediments  before 
and  after  removing,  by  chemical  attack,  vari- 
ous components  from  the  sediments.  The  ma- 
jor decrease  in  sorption  was  observed  after  re- 
moval of  organic  matter  (extraction  with 
NaOH);  they  thus  attributed  an  important  role 
to  organic  matter  in  the  sorption  of  Cd. 

In  line  with  the  above  findings,  we  assume 
here  that  the  two  main  components  of  oxic 
sediments  that  sorb  Cd  are  Fe  oxyhydroxides 
and  organic  matter,  that  these  components 
have  surface  sites  for  sorption  that  can  be  treat- 
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ed  as  surface  ligands,  and  that  surface  com- 
plexation  concepts  developed  for  simple,  well- 
defined  systems  in  the  laboratory  can  be 
applied  to  more  complex  sediment  compo- 
nents. According  to  surface  complexation  the- 
ory, the  adsorption  of  Cd  on  Fe  oxyhydroxide 
surfaces  can  be  expressed  (Benjamin  and  Leck- 
ie  1981)  by 


=  Fe-OH,  +  Cd=^ 


'K, 


=  Fe-OCd-l-;cH+  (2) 


=  Fe-OCd}[H"]^ 
=  Fe-OH,}[Cd-*] 


(3) 


where  charges  on  the  solid  species  are  omitted 
for  simplicity,  */^F<ï-cd  is  an  apparent  overall 
equilibrium  constant,  x  is  the  average  apparent 
number  of  protons  released  per  Cd^*  ion  ad- 
sorbed (Honeyman  and  Leckie  1986),  and 
{^Fe-OH,}  and  (  =  Fe-OCd}  represent  re- 
spectively the  concentration  of  free  surface  sites 
on  the  Fe  oxyhydroxides  and  the  concentra- 
tion of  sites  occupied  by  Cd.  Throughout  this 
paper,  the  notation  "  =  "  refers  to  adsorption 
sites,  whereas  { • }  and  [  •  ]  refer  to  concentra- 
tions of  solid  and  dissolved  species  respec- 
tively. 

At  low  adsorption  density  (i.e.  when  the  con- 
centration of  occupied  adsorption  sites  is  small 
compared  to  the  free  site  concentration)  the 
condition 


=  Fe-0-}T  =  {  =  Fe-OHJ 


(4) 


should  apply,  where  {^Fe-O-lx  is  the  total 
concentration  of  sites.  This  latter  concentra- 
tion can  in  turn  be  expressed  as 

{  =  Fe-0-}T  =  iVFc  X  {Fe-ox}  (5) 

where  Nf^  is  the  number  of  moles  of  adsorption 
sites  of  the  Fe  oxyhydroxides  per  mole  of  Fe, 
and  {Fe-ox}  is  the  analytical  concentration  of 
Fe  oxyhydroxides.  If  it  is  assumed  that  only 
one  site  is  occupied  by  each  sorbed  metal  ion, 
then 

(=Fe-OCd}  =  {Fe-Cd)  (6) 

where  {Fe-Cd}  is  the  analytical  concentration 
of  Cd  associated  with  the  Fe  oxyhydroxides. 
Combining  Eq.  3  to  6  leads  to 


{Fe-Cd} 


where  AV^-cd  is  an  apparent  equilibrium  con- 
stant which  is  a  function  of  pH. 

By  similar  reasoning,  the  following  expres- 
sion can  be  derived  for  sorption  of  Cd  on  or- 
ganic matter  (OM): 


K^. 


*AoM- 


[H-]" 

{OM-Cd} 

{OM}[Cd=*] 


(8) 


{ Fe-ox  }[Cd2-^] 


(V) 


where  by  analogy  with  Eq.  2-7,  N^^  is  the  num- 
ber of  moles  of  sites  on  the  organic  matter  per 
mole  organic  C,  and  Aro^,_c<j  and  *AoM-cd  are 
apparent  overall  equilibrium  constants.  Sim- 
ilarly, y  is  the  apparent  number  of  protons 
released  per  Cd^*  ion  adsorbed  on  organic 
matter,  and  {OM}  and  {OM-Cd}  are  the  con- 
centrations of  organic  matter  and  of  Cd  as- 
sociated with  organic  matter  respectively.  Lin- 
earization of  Eq.  7  and  8  yields 

log  ^Fe^a  =  xpH  +  \og(N,,  X  «Af,^,)      (9) 

log  A:„^d  =  y  pH  -I-  log(jVoM  X  *Ko^i)       (  1 0) 

and  plots  of  log  ATfc-cci  or  log  K^^^a  vs.  pH 
should  yield  slopes  of  x  and  y  and  intercepts 
on  the  y-axis  of  log(A^Fe  ^  *'^Fe-cd)  and  log(A^oM 
^  */Q,„-cd>  respectively. 

We  have  estimated  Kr^_ç^  and  K^^^-ca  for  each 
sampling  site  from  determinations  of  the  three 
variables  on  the  right-hand  side  of  Eq.  7  and 
8.  In  the  calculation  of  the  first  of  these  vari- 
ables, [Cd^"^]  (nmol  liter"';  Table  2),  only  in- 
organic ligands  were  considered;  because  of  the 
lack  of  relevant  equilibrium  constants,  possi- 
ble Cd  complexation  by  natural  organic  li- 
gands could  not  be  considered.  The  calculated 
value  of  [Cd-*]  in  these  waters  is  very  close  to 
that  of  [Cd];  calculated  contributions  of 
CdS04(aq)  and  CdCr  are  generally  <5%,  the 
only  exceptions  being  for  lakes  St.  Charles 
(10%),  McFarlane  (15%)  and  Silver  (19%),  all 
in  the  Sudbury  area. 

Estimation  of  the  concentrations  of  Cd  as- 
sociated with  the  Fe  oxyhydroxides  or  with 
organic  matter  ({ Fe-Cd  1  and  {OM-Cd}),  two 
of  the  remaining  variables  in  Eq.  7  and  8,  is 
however  not  straightfonAard.  Most  of  sedi- 
ment Cd  was  found  in  extracts  (1)  35 ±20%,  ' 
(2)  29 ±10%,  (3)  15±10%,  and  (4)  23 ±20%; 
concentrations  of  Cd  in  extract  (5)  were  de- 
tectable only  in  a  fev.  cases  and  levels  in  extract 
(6)  were  undetectable  for  all  sediment  samples. 
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Large  proponions  of  sedimeniary  Cd  were  thus 
extracted  with  relatively  mild  reagents.  In  such 
a  situation,  it  is  difficult  to  decide  objectively 
which  portion  of  this  easily  extractable  Cd  was 
associated  with  Fe  oxyhydroxides  or  with  or- 
ganic matter  in  the  onginal  sediment. 

To  overcome  this  difficulty,  we  proceeded 
in  the  following  manner.  The  ratio  \Cd]^: 
{Felc,  obtained  from  dissolution  of  the  Fe  ox- 
ide deposit  on  the  Teflon  collectors,  was  as- 
sumed to  be  representative  of  the  ratio  {Fe- 
Cd}  :  {Fe-ox}  prevailing  in  the  sediments.  The 
values  of  [Cd-"],  together  with  those  of  the 
ratios  {Cd}^:  {Fel^  obtained  from  the  Teflon 
collectors,  were  then  used  to  calculate  /Cpc-cd 
according  to  Eq.  7.  The  plot,  corresponding  to 
Eq.  9,  is  shown  in  Fig.  2 A;  linear  regression  of 
these  data  yields 

log  ATpc-cd  =  0.82(±0.08)  pH  -  l.30(±0.38) 
(n=  16;         r=  =  0.89).  (H) 

The  total  sorbed  Cd  concentration  is  written 
{Cd},  =  (Fe-Cd}  +  (OM-Cd}       (12) 

since  we  assume  that  only  Fe  oxyhydroxides 
and  organic  matter  bind  Cd  in  sediments. 
Combining  Eq.  7  and  12  allows  calculation  of 
{OM-Cd}  according  to 


6  - 


{OM-Cd} 
=  {Cd},- 


Nf.  x*/CFc^d{Fe^x}[Cd^-]     ^j3^ 


[H-]' 


In  the  calculation  of  {OM-Cd}  with  Eq.  1 3, 
the  value  of  {Fe-ox}  (molFeg"'  sediment  dry 
wt;  Table  2)  was  taken  as  the  sum  of  the  Fe 
concentrations  extracted  from  the  sediments 
in  steps  (3)  and  (4)  with  the  reducing  reagent 
NH.OH  •  HQ.  Measurements  with  pure  phases 
prepared  in  the  laboratory  have  shown  that 
this  reagent  dissolves  synthetic  ferrihydrite  and 
lepidocrocite  almost  quantitatively  (Fortin  et 
al.  in  prep.),  i.e.  the  Fe  oxyhydroxide  forms 
that  should  be  most  effective  in  binding  Cd. 
Together,  these  two  extractions  dissolved  most 
(89  ±7%)  of  the  Fe  extractable  in  the  first  five 
extractions.  Values  of. X  (0.82)  and  AVe  ^  *^F<s-cd 
(10-"°)  were  taken  from  Eq.  II.  Values  of 
{OM-Cd}  calculated  on  this  basis  with  Eq.  13 
were  used,  together  with  [Cd=-]  and  {OM} 
(Table  2),  to  calculate  K^^^^  according  to  Eq. 
8.  The  values  of  sedimentary  organic  C  con- 
centrations (mol  organic  C  g"'  sediment  dry 
wi)  were  used  as  estimates  of  {OMJ.  The  plot 


Fig  2  Apparent  overall  equilibrium  constants  for  the 
sorption  of  Cd  on  natural  Fe  oxyhydroxides  and  organic 
mailer  as  a  function  of  pH.  Symbols  as  in  Fig.  1. 


corresponding  to  Eq.  10  is  shown  in  Fig.  2B 
and  regression  of  the  data  yields 
logJ^oM^d  =  0.97(±0.09)  pH  -  2.45(±0.46) 
(^j  =  35;        r^  =  0.78).  (M) 

Sites  located  in  lakes  Gullfeather  and  Bigwind, 
for  which  the  pore-water  profiles  showed  ev- 
idence of  anoxic  conditions  at  the  sediment- 
water  interface  (release  of  dissolved  Fe  from 
the  sediment;  reduction  of  sulfate  at  the  inter- 
face), were  not  included  in  the  calculation;  sta- 
tion MA-01,  consistently  an  outlier,  was  also 
excluded. 

LogA:F^.  and  logATo^-cd  show  pH  depen- 
dency (Fig.  2),  as  expected  from  Eq.  9  and  10; 
this  behavior  is  consistent  with  the  simple  sur- 
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face  complexation  model  used,  despite  the 
crude  assumptions  that  were  made.  Many  fac- 
tors may  contribute  to  the  scatter  evident  in 
Fig.  2.  The  simple  model  retained  does  not 
consider  electrostatic  interactions  for  the  sorp- 
tion of  Cd  on  Fe  oxyhydroxides  and  organic 
matter,  i.e.  it  assumes  that  the  free  energy  of 
adsorption  is  dominated  by  chemical  inter- 
actions rather  than  by  electrostatic  interac- 
tions. Values  of  the  parameters  A^pc.  j'^om. 
*A'Fj-cd>  and  *Ko^,_(~ci  are  considered  here  to  be 
constants,  whereas  they  might  be  expected  to 
vary  with  the  nature  of  sedimentary  Fe  oxy- 
hydroxides and  organic  matter  (Luoma  and 
Davis  1983);  it  is  probable  that  the  nature  of 
these  two  components  present  in  the  top  0.5 
cm  of  the  lake  sediments  varies  among  sta- 
tions. Extraction  of  Fe  with  NH^OHHCl  is 
Only  moderately  selective  and  will  solubilize 
various  Fe  oxyhydroxide  forms;  similarly,  to- 
tal organic  sedimentary  C  is  a  very  crude  es- 
timator of  the  organic  matter  active  in  the 
sorption  process.  Only  two  sediment  compo- 
nents, Fe  oxyhydroxides  and  organic  matter, 
have  been  assumed  to  be  responsible  for  sorp- 
tion; sorption  by  Mn  oxyhydroxides,  which 
might  be  important  in  some  high  pH  lakes,  has 
been  neglected.  According  to  laboratory  ex- 
periments in  well-defined  systems,  the  sorp- 
tion constants  should  also  vary  with  the  den- 
sity of  adsorption  (Benjamin  and  Leckie  198 1) 
and  with  the  concentration  of  particles  (Di- 
Toro  et  al.  1986),  both  of  which  vary  among 
the  sites.  Sedimentary  Co  concentrations 
probably  do  not  vary  as  rapidly  as  dissolved 
Cd  and  H"^  in  response  to  variations  in  envi- 
ronmental conditions,  i.e.  the  assumed  sorp- 
tion equilibrium  might  not  be  fully  attained. 
The  concentration  (and  possibly  the  nature)  of 
the  dissolved  organic  matter  varies  among  sta- 
tions (Table  2);  in  calculating  Kp^^^  and  K>^*-cd, 
possible  complexation  of  Cd  by  dissolved  or- 
ganic matter  was  not  taken  into  account. 

Some  of  the  scatter  in  Fig.  2  might  also  be 
due  to  experimental  difficulties.  The  dissolved 
Cd  concentrations  are  low,  especially  for  high 
pH  lakes;  such  samples  are  particularly  subject 
to  contamination  and  the  analytical  determi- 
nations themselves  are  less  precise  than  at 
higher  concentrations.  During  collection  of  the 
surficial  sediment,  the  oxic  layer  may  be  con- 
taminated by  reduced  sediments,  leading  to 
incorrect  estimation  of  {Cd}r,  {Fe-ox},  and 
(OM|  present  in  the  surficial  sediment.  The 


ratio  (Cd}^:  {Fe}^  is  also  assumed  to  be  rep- 
resentative of  the  ratio  {Fe-Cd}  :  {Fe-ox}  pre- 
vailing in  the  oxic  sediments.  The  Teflon  sheets, 
however,  collect  "young"  diagenetic  Fe  oxy- 
hydroxides; for  example,  no  goethite  could  be 
identified  in  the  Fe-rich  material  collected 
(Fortin  et  al.  in  prep.),  whereas  there  are  in- 
dications that  this  crystalUne  solid  can  be  pres- 
ent in  lake  sediments  (Schwertmann  et  al. 
1987).  In  addition,  chemical  analysis  of  the 
Fe-rich  diagenetic  deposits  indicated  the  pres- 
ence of  organic  C  and  Mn  (Fortin  et  al.  in 
prep.);  organic  matter  and  Mn  oxyhydroxides 
(for  higher  pH  lakes)  could  be  responsible  for 
binding  a  fraction  of  the  Cd  collected  on  the 
Teflon  sheets  and  assumed  to  be  bound  to  Fe 
oxyhydroxides. 

Comparison  of  {Cdjr  and  {OM-Cd}  in  Ta- 
ble 2  indicates  that  the  association  of  Cd  in 
oxic  lake  sediments  is  dominated  by  its  inter- 
actions with  organic  matter.  For  the  lake  sed- 
iments studied,  {OM-Cd}  and  {Fe-Cd} 
represent  92±8%  and  8±8%  of  {Cd}x,  re- 
spectively. The  proportion  of  each  form  varies 
from  one  site  to  another,  depending  on  the 
relative  concentrations  of  sedimentary  Fe  oxy- 
hydroxides and  organic  matter  in  the  sediment 
and  on  lake  pH.  It  sholild  be  noted  that  similar 
conclusions  were  reached  by  Luoma  (1986)  for 
Cu  in  oxidized  estuarine  sediments.  The  val- 
ues of  {OM-Cd}  calculated  with  Eq.  13  cor- 
respond very  closely  to  the  sum  of  Cd  concen- 
trations found  in  sediment  extracts  (1),  (2),  (4), 
and  (5)  (slope  =  0.95;  r^  =  0.99;  n  =  35).  Cal- 
culated values  of  {Fe-Cd}  agree  also  reason- 
ably well  with  Cd  concentrations  found  in  the 
remaining  extract  (3):  slope  =  1.26;  r^  =  0.65; 
n  =  35.  It  is  however  difficult,  at  the  moment, 
to  rationalize  these  correlations  on  a  geochem- 
ical  basis. 

The  model  developed  above  should  prove 
useful  (if  our  assumptions  regarding  the  cal- 
culation of  A^oM-cd  are  valid)  for  obtaining  a 
rough  estimate  of  the  partitioning  of  Cd  be- 
tween organic  matter  and  Fe  oxyhydroxides  in 
oxic  sediments  of  other  lakes.  Combination  of 
Eq.  7,  8,  and  12  leads  to 

{Fe-Cd} 


{Cd}, 
=  (AV.  X  '/^p,^, { Fe-ox  [H^]') 

-^(Np,  X  •/:p,^,.:Fe^xf[H-]' 

+  A',.,  X  ''K.,^,^^.OM  [H'Y)  (15) 
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(OM-Cdl 


(Cd!, 


^  •A;..^jiOM}[H*]') 


+  A'o„  X  •A;,^,{0M}[H*]').  (16) 

The  only  requirement  would  be  to  measure 
{Cdlr  (sum  of  nondetrital  forms),  {OM},  and 
(Fe-ox}  for  the  same  oxic  sediment  sample 
and  lake  pH;  the  necessary  values  of  x,  y,  N^^ 
X  *AV,^o.  and  A'o„  x  *Ao„_cd  are  obtained 
from  linear  regressions  1 1  and  14.  The  values 
of  iOM-Cd}  calculated  with  Eq.  16  for  the 
lakes  studied  are  shown  in  Table  2. 

Relationship  between  Cd  concentrations  in 
water  and  in  A.  grandis  — Cd  concentrations  in 
various  soft  tissues  of  the  freshwater  bivalve 
A.  grandis  are  given  in  Table  3.  Concentrations 
in  the  whole  organism,  [Cd(Org)],  were  cal- 
culated as  follows: 


[Cd(Org)]  = 


2[Cd(tissue)],  W, 


(17) 


Tabic  3.  Cd  concentrations  (>ig  g  ')  in  tht  tissues  of 
Anodonia  grandis  from  vanous  sampling  sites.  He  — hc- 
patopancreas;  Gi  — gills;  Ma  — mantle:  Re  — remaining; 
Org— reconstituted  organisms  (calculated  with  Eq.  17). 


where  [Cd(tissue)],  and  W^,  are  the  Cd  concen- 
tration and  the  dry  weight  of  the  /th  tissue. 
Table  3  shows  a  50-fold  variation  in  [Cd(Org)] 
among  stations.  Concentrations  are  generally 
highest  in  the  gills  of  the  bivalve,  especially 
for  the  more  polluted  sites.  Gills  always  con- 
tained the  greatest  proportion  (40  ±  1 3%)  of  the 
total  Cd  burden  of  whole  animals;  the  mantle, 
hepatopancreas,  and  remaining  tissues  con- 
tributed 21  ±7,  ll±8,and28±l  1%.  Stepwise 
multiple  regression  was  performed  to  examine 
the  relationships  between  [Cd(Org)]  as  the  de- 
pendent vanable  and  [Cd'*],  [Ca-*],  [H*],and 
Cors  as  the  independent  variables.  Due  to  in- 
tercorrelations  among  both  pH  and  Co^g  and 
pH  and  [Cd-*],  the  largest  partial  regression 
coefficient,  found  for  [Cd-*],  explained  most 
of  the  variance  (r-  =  0.80);  addition  of  pH  (P 
<  0.08)  increased  the  explained  variance  by 
only  3%  (H  =  0.83). 

Figure  3  illustrates  the  dependence  of 
[Cd(Org))  on  the  concentration  of  the  free  Cd-  * 
ion,  as  calculated  from  [Cd]  and  inorganic  li- 
gand  concentrations.  The  linear  regression 
equation  found  is 

[Cd(Org)]  =  44(  =  5)  [Cd=*]  -^  10(±15) 

{n  =  17;         H  =  0.81)  (18) 

with  [Cd-*]  in  nmol  liter"'.  Regressions  be- 
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tween  [Cd^*]  and  Cd  concentrations  in  the  tis- 
sues are  highly  significant  (P  <  0.001),  except 
for  Cd  in  the  hepatopancreas  which  is  signif- 
icant (Table  4). 

Laboratory  bioassays,  conducted  with  ma- 
rine benthic  animals  exposed  to  dissolved  met- 
als under  carefully  controlled  conditions,  have 
shown  that  the  effects  of  trace  metals  are  re- 
lated to  the  free-metal  ion  concentration  and 
not  to  the  total  metal  concentrations.  For  ex- 
ample, mortality  of  the  shrimp  Palaemonetes 
pugio  exposed  to  Cd  was  shown  to  be  propor- 
tional to  [Cd^*]  (Sunda  et  al.  1978).  Similarly, 
shon-term  (14  d)  accumulation  of  Cu  in  the 
oyster  Crassostrea  virginica  was  found  to  be 
proportional  to  [Cu-*]  (Zamuda  and  Sunda 
1982).  These  laboratory  observations  are  con- 
sistent with  our  field  results  (Fig.  3;  Table  4). 

However,  we  would  have  obtained  similar 
regression  equations  with  [Cd]  since  it  is  close 
to  the  computed  value  for  [Cd^*]  at  all  sites. 
Given  the  similarity  between  [Cd-*]  and  [Cd] 
values,  and  the  relative  uncertainty  in  the  re- 
lations given  in  Table  4.  we  cannot  demon- 
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Fig.  3.  Relationship  between  Cd  concentrations  in  (he 
tissues  of  Anodonta  grandis  and  [Cd-*]  calculated  from 
total  dissolved  Cd  and  inorganic  ligands.  Symbols  as  in 
Fig.  1. 

Strate  unambiguously  that  the  free-ion  activity 
model  (Morel  1983)  applies  to  Cd  accumula- 
tion in  A.  grandis.  This  model  has  been  shown 
to  apply  to  various  organisms  and  various 
metals  under  laboratory  conditions  where  ar- 
tificial ligands  were  used,  but  has  not  been  val- 
idated in  the  field.  It  is  unlikely  that  A.  grandis 
will  be  found  in  freshwaters  where  inorganic 
ligands  (SO4,  CI)  are  present  in  sufficiently  high 
concentrations  to  complex  Cd  significantly;  an 


alternative  approach  10  demonstrate  applica- 
tion of  the  model  to  Cd  accumulation  by  A. 
grandis  in  the  field  would  be  lo  find  sites  where 
natural  organic  matter  complexes  Cd  signifi- 
cantly, provided  we  can  account  for  this  com- 
plexation  and  correctly  estimate  [Cd'*].  Par- 
tial correlation  analysis  (holding  Cd  constant) 
did  not  show  any  significant  relationship  be- 
tween [Cd(Org)]  and  DOC  (r  =  -0. 1 3;  «  =^  1 7), 
suggesting  that  complexaiion  by  organic  mat- 
ter did  not  influence  Cd  uptake  by  the  bivalves 
at  our  sites. 

Reasonable  predictions  of  Cd  concentra- 
tions in  the  tissues  of  .'1.  grandis  can  be  ob- 
tained presently  from  ambient  [Cd-*]  values 
with  the  regression  equations  given  in  Table 
4.  The  predictive  power  of  these  equations 
could  however  be  improved.  Much  of  the  scat- 
ter in  Fig.  3  is  probably  due  to  changes  in 
dissolved  Cd  on  short  time  scales  with  envi- 
ronmental conditions  (phytopiankton  bio- 
mass,  pH,  etc.).  For  example,  Yan  et  al.  (  1 990) 
observed  a  twofold  change  in  [Cd]  in  Red  Chalk 
Lake  during  the  ice-free  season;  this  lake  is 
within  our  study  area  (Muskoka).  The  bivalves 
would  not  react  rapidly  to  such  changes  in  [Cd]. 
Indeed,  we  have  observed  that  Cd  concentra- 
tions in  the  tissues  of  A.  grandis  specimens 
transplanted  from  unpolluted  Lake  Brompton 
to  polluted  Lake  Joannes  (both  included  in  our 
present  study)  had  not  reached,  after  3  yr,  the 
Cd  levels  observed  in  the  indigenous  speci- 
mens of  Lake  Joannes  (Tessier  unpubl.  re- 


Table  4.  Linear  regression  equations  describing  Cd  accumulation  in  various  tissues  of  A  nodonla  grandis  as  a  function 
of  [Cd]  or  sedimentary  variables.  Tissues  abbreviations:  Gi— gills;  Ma— mantle;  He  — hepatopancreas;  Re— remains; 
Org— whole  organism.  ('*•:  P  <  0.001;  **:  P  <  0.01.) 


Dissolved  Cd 
[Cd(tissue))  =  F  [Cd-'*]  +  [CdCtissue)]" 
[Cd(Gi)]=  l08(±l5)[Cd'*]  +  I6(±41) 
[Cd(Ma)]  =  36(±6)  [Cd^*]  +  17(±16) 
(Cd(He)]  =  27{±8)  [Cd=*]  +  16(±21) 
[Cd(Re)]  =  21(±3)  [Cd=*l  +  5(±9) 
[Cd(Org)l  =  44(±5)[Cd-*]+  10(±15) 

Sedimentary  vanables 

,^     .         ,  ^ FiCdU[H']'"' 

I     (  issue^j       ,o-i>o,pe_ox)[H*]""  +  lO-^-'MOMKH- 

[Cd(Gi)l  =  123(±24)  [Cd=*]  +  28(±64) 

[Cd(Ma)]  =  52(±5)  [Cd-*]  +  14(±14) 

[Cd(He)]  =  41(±6)  [Cd-*]  -i-  13(±15) 

[Cd(Re)]  =  33(±3)  [Cd=*l  -i-  5(±8) 

[Cd(Org)]  =  59(±7)  [Cd-*]  +  I  l(±  18) 


«=17 

r-  =  0.77 

/!=  17 

r"  =  0.71 

/!=    17 

r-  =  0.45 

n=  17 

r-  =  0.74 

n  =  17 

r  =0.81 

[Cd(tissue)]" 

n  =  19;  r-  =0.61 

/!=  19 

r  =  0.86 

/7=  19 

r'  =  0.77 

n=  19 

r-  =  0.87 

n=  19 

r-  =  0.82 
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suits).  The  use  of  iimc-avcragcd  values  of 
[Cd-*]  (or  [Cd])  in  the  linear  regression  equa- 
tions, instead  of  the  single  values  used  here, 
would  undoubtedly  improve  the  predictive 
power  of  these  equations  for  Cd  concentrations 
in  the  tissues  of  .-1.  grandis.  Measurements  of 
temporal  changes  m  Cd  concentrations  in  the 
tissues  of  transplanted  bivalves  suggest  that 
they  accumulate  metals  only  dunng  the  warm 
season  (e.g.  from  May  to  October).  For  opti- 
mum predictive  power,  the  mean  [Cd-*]  val- 
ues should  then  be  obtained  preferentially  for 
this  penod.  Another  source  of  scatter,  proba- 
bly less  important,  is  variability  in  Cd  con- 
centration in  the  tissues  of  the  bivalves  due  to 
biological  factors  (length  or  age,  reproductive 
condition);  in  the  present  study,  we  attempted 
to  minimize  the  influence  of  such  factors  by 
collecting  specimens  in  a  narrow  size  range,  at 
times  chosen  to  avoid  the  presence  of  gravid 
animals. 

Relationship  between  Cd  in  the  sediments 
and  in  A.  grandis  — Since  [Cd-*]  is  difficult  to 
measure,  an  alternative  approach  to  the  use  of 
dissolved  Cd  to  predict  [Cd(Org)l  would  be  to 
consider  sedimentary  Cd  concentrations  as 
possible  predictors.  When  Cd  concentrations 
in  the  bivalves  (Table  3)  are  compared  with 
{Cd},  in  the  sediments  (Table  2),  no  relation- 
ship is  shown  (Fig.  4).  The  same  conclusion  is 
reached  when  the  Cd  concentration  in  the 
whole  organism  or  any  soft  tissue  is  compared 
with  Cd  concentrations  found  in  any  sediment 
extract,  whether  or  not  the  extraclable  Cd  is 
normalized  with  respect  to  sediment  Fe  oxide 
or  organic  C  concentrations.  Such  normal- 
izations, based  on  sound  geochemical  princi- 
ples regarding  the  binding  strength  of  sediment 
components  toward  trace  metals  (Luoma  and 
Bryan  1978),  have  been  useful  for  predicting 
As,  Cu,  Hg.  and  Pb  bioaccumulation  in  bi- 
valves (Langston  1980,  1982;  Luoma  and  Bry- 
an 1978;Tessieretal.  1983.  1984),  but  do  not 
seem  adequate  for  predicting  Cd  accumulation 
in  A.  grandis  from  a  large  variety  of  lakes  as 
in  the  present  study. 

Close  examination  of  the  data  in  Tables  2 
and  3  reveals  that,  at  a  given  ratio  of  {Fe-Cd}  : 
{Fe-ox:  or  'OM-Cd}  :  {OM},  the  bivalves  ac- 
cumulate much  more  Cd  in  their  tissue  in  lakes 
of  low  pH  than  in  those  of  high  pH.  For  ex- 
ample, at  site  VA-01  ({Fe-Cd;  :  i Fe-ox}  =  2.7 
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Fig.  4.  Relationship  between  Cd  concentrations  in  the 
tissues  oi  Anodoma  grandis  and  (Cdl,  in  the  sediments. 
Symbols  as  in  Fig.  1. 


X  10"*;  pH  =  6.57)  the  bivalves  accumulated 
152  Mg  g~'  of  Cd  whereas  at  site  BR-05  ({Fe- 
Cd}  :  (Fe-ox}  =  3.4  x  lO"*;  pH  =  7.33)  they 
accumulated  only  40  Mg  g~  '  ■  Similarly,  the  Cd 
concentration  in  the  bivalves  was  much  higher 
atsiteBO-02((OM-Cd}:{OM}  =8.8  x  lO-*; 
pH  =  6.34;  [Cd(Org)]  =  1 10  ^g  g"')  than  at 
site  H-01({OM-Cd}:(OM}  =  8.1  x  10"*;  pH 
=  6.93;  [Cd(Org)]  =  15  Mg  g"').  This  pH  de- 
pendence of  Cd  bioaccumulation  at  a  given 
concentration  ratio  of  sorbed  Cd  :  sorbent  can 
be  taken  into  account  by  combining  the  surface 
complexation  concepts  described  above  with 
the  free-metal  ion  activity  model,  as  described 
below. 

If  the  Cd  concentration  in  the  molluscs  is 
assumed  to  be  related  to  [Cd^*],  as  suggested 
by  Fig.  3,  one  can  write 

[Cd(Org)]  =  F[Cd=-]  +  [Cd(Org)]o,  (19) 

which  is  of  similar  form  to  the  linear  regression 
18;  /"  is  a  proportionality  factor  between 
[Cd(Org)]  and  [Cd^*]  whereas  [Cd(Org)]°  is  the 
intercept  on  the  >-axis.  Combining  Eq.  7,  8, 
12,  and  19  vields 


[Cd(Org)]  =  F- 


(  Fe-Cd;  [H-]' 


AV.  X  •A:f,_c<j{ Fe-ox} 
[Cd(Org)]^ 


(20a) 
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Fig.  5.  Relationship  between  Cd  concentrations  in  the 
\\%%\it^oî Anodoma grandis zr^A  the  [Cd^*]  calculated  with 
Eq.  21c.  Symbols  as  in  Fig.  I. 


[Cd(Org)]  =  F 


{OM-Cd}[H-^]>- 
+  [Cd(Org)]o, 


(20b) 


and 
[Cd(Org)]  =  F 


(Cd},[H+]^+> 


+  [Cd(Org)]°. 


(20c) 


The  right-hand  quotients  in  Eq.  20a,b,c  can 
be  used  to  evaluate  [Cd^"^]  without  need  to 
measure  [Cd],  provided  that  the  geochemical 
constants  X.  y.  A^F,  x  •ATp^.^d,  and  N^^  x  *K^„,.j^ 
are  known.  All  three  quotients  lead  to  the  same 
value  of  [Cd^*].  Substituting  the  values  of  these 
constants  obtained  from  linear  regressions  1 1 
and  14  into  Eq.  20,  one  obtains 


[Cd(Org)]  =  F 


[Fe-Cd}[H^]°" 
10-'^°{Fe-ox} 


+  [Cd(Org)]o, 


[Cd(Org)]  =  F 


and 


{OM-Cd}[H-^f^^ 
■   lO-^-'MOM} 
[Cd(Org)]o, 


(21a) 


(21b) 


[Cd(Org)]  =  F 


(Cd},[H*]' 


(10-'5°(Fe-ox}[H-^]0" 
+  IQ-^-'MOMKH^lo'^), 

+  [Cd(Org)]o.  (21c) 

Similar  equations  can  be  written  for  indi- 
vidual tissues  of  the  organisms.  Equation  21 
can  be  used  in  practice  to  predict  Cd  concen- 
trations in  the  tissues  of  various  aquatic  or- 
ganisms. The  first  step  is  to  determine  the  val- 
ues of  F  and  [Cd(Org)]°.  We  have  determined 
these  two  constants  with  Eq.  21c  for  various 
tissues  of  ^.  grandis  with  the  values  of  {Cd}T, 
pH,  {Fe-ox},  and  {OM}  given  in  Table  2  and 
the  Cd  tissue  concentrations  given  in  Table  3 
{see  Table  4;  Fig.  5).  We  found  that  using  the 
sum  of  Cd  extracted  in  steps  (1),  (2),  (3),  and 
(4)  instead  of  {OM-Cd}  in  Eq.  21b  did  not 
change  the  explained  variance  (H).  Substitu- 
tion of  iOM-Cd}  by  {Cd}^  in  Eq.  21b  reduced 
the  value  of  r^,  but  only  slightly  (e.g.  from  0.82 
to  0.79  for  [Cd(Org)]).  This  lack  of  sensitivity 
reflects  the  fact  that  the  binding  of  Cd  to  the 
sediment  is  dominated  by  organic  matter. 

The  values  of  Fand  [Cd(Org)]°  obtained  from 
the  linear  regressions  could  be  used  to  predict 
Cd  concentrations  in  the  tissues  of  A.  grandis 
in  other  lakes,  provided  that  lake  pH  and  {Cdlr, 
{Fe-ox},  and  {OM}  have  been  determined  for 
the  same  sediment  samples  from  each  lake. 
Equation  2 1  may  be  valid  to  predict  Cd  con- 
centrations in  organisms  other  than  A.  grandis 
that  obtain  Cd  from  the  water  (or  from  food 
such  as  phytoplankton,  the  Cd  concentration 
of  which  would  be  related  to  that  in  the  water); 
the  only  requirement  would  be  to  determine 
the  values  of  F  and  [Cd(Org)]°,  which  should 
be  organism-specific.  Because  the  geochemical 
constants -v.y,  Np,  x  «Ap,^,  and  A^o„  x  *K„„_^ 
are  involved  in  predicting  Cd  accumulation  by 
various  organisms,  their  precise  determination 
would  yield  a  better  basis  for  predicting  Cd 
bioavailability  to  benthic  organisms. 

Use  of  sedimentary  measurements  as  in  Eq. 
21  to  predict  Cd  bioaccumulation  offers  con- 
siderable advantages  over  direct  measurement 
of  dissolved  Cd  (as  in  Eq.  18).  Measuring  low 
values  of  [Cd]  in  freshwaters  is  difficult  and 
requires  the  use  of  trace  metal-free  techniques; 
possibilities  of  inadvertent  contamination  are 
rampant.  An  additional  problem  is  the  tem- 
poral vanation  of  dissolved  Cd;  considerable 
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cirort  would  be  needed  to  obtain  representative 
mean  values  of  total  dissolved  Cd  at  each  sam- 
pling site.  By  companson,  measunng  iCd},, 
|Fe-ox},  {OMI,  and  lake  pH  is  simpler. 

In  earlier  studies,  field  measurements  along 
trace  metal  gradients  have  shown  that  predic- 
tion of  metal  concentrations  in  estuarine  and 
freshwater  bivalves.  [M(Org)].  was  greatly  im- 
proved when  the  trace  metal  concentration  ex- 
tracted from  the  sediments  was  normalized 
with  respect  to  the  Fe  oxide  or  organic  matter 
content  of  the  sediments.  Hence,  the  ratios 
{Pb}  :  (Fe}  and  {As}  :  {Fe},  all  extracted  with 
1  N  HCl.  and  the  ratio  {Hg  extracted  with 
HNOj}  :  {OM\  were  found  to  be  the  best  pre- 
dictors of  Pb,  As,  and  Hg  concentrations  re- 
spectively in  the  estuarine  bivalve  Scrobicu- 
laria  plana  {Langsion  1980,  1982;  Luoma  and 
Bryan  1978).  Similarly,  the  best  predictors  of 
Cu  and  Pb  in  the  tissues  of  the  freshwater  bi- 
valves A.  grandis  and  Elliplio  complanala  were 
found  to  be  the  ratios  {Cu}  :  {Fef  and  {Pb}  : 
{Fe}  extracted  with  NH.OHHCl,  a  reducing 
reagent  (Tessier  el  al.  1983,  1984).  In  other 
words,  according  to  these  studies,  trace  metal 
concentrations  in  the  organisms  are  best  ex- 
pressed as 


[M(Org)]  =  k 


{S„-M] 


[M(Org)l°     (22) 


where  A:  is  a  proportionality  constant,  {5,}  the 
concentration  of  sediment  component  n  (which 
can  be  Fe  oxyhydroxides  or  organic  matter), 
{5„-M}  the  concentration  of  the  trace  metal 
M  associated  with  sediment  component  n.  and 
(M(Org)]°  the  intercept  on  the  v-axis  (usually 
close  to  zero). 

The  present  study  offers  a  likely  explanation 
for  these  findings.  The  studies  of  Langston 
(1980,  1982)  and  Luoma  and  Bryan  (1978) 
were  performed  in  estuaries,  where  the  pH  is 
relatively  constant;  the  freshwater  studies 
(Tessier  et  al.  1983,  1984)  were  earned  out  in 
three  lakes  located  in  a  restncted  geographical 
area  where  the  lake  pH  is  relatively  constant. 
In  such  cases,  [H*]'  and  [H']'  become  ap- 
proximately constant  and  Eq.  20a  and  20b  re- 
duce to  Eq.  22.  The  success  of  sediment  Fe  or 
organic  C  as  the  normalizing  factor  in  these 
studies  could  be  due  to  its  dominance  in  sorp- 
tion of  the  particular  trace  metal  investigated. 
For  example,  the  strong  association  of  Pb  (Lion 


ctal.  1982;  Balistnen  and  Murray  1982;  Leck- 
le  et  al.  1984)  and  As  (Aggett  and  Roberts 
1 986;  Belzile  and  Tessier  1 990)  with  sedimen- 
tary- Fe  oxyhydroxides  is  well  documented; 
similarly,  the  affinity  of  Hg  for  sediment  or- 
ganic matter  (NRCC  1 988)  is  well  known.  Thus 
the  simpler  normalization  approach  appears 
to  be  applicable  only  where  a  pH  gradient  does 
not  exist. 

From  the  preceding  discussion,  it  follows 
that  the  observation  of  a  strong  relationship 
between  trace  metal  levels  in  a  benthic  species 
and  sediment  charactenstics  (e.g.  extraciable 
metal  concentration  normalized  with  respect 
to  the  concentration  of  a  sediment  component 
like  Fe  oxyhydroxides  or  organic  matter)  can- 
not be  construed  as  evidence  that  the  main 
route  of  trace  metal  uptake  is  via  ingestion  of 
sediment  particles.  Such  relationships  could 
also  be  observed  for  organisms  that  obtain 
metal  directly  from  the  water  (or  from  food, 
the  Cd  content  of  which  would  be  related  to 
that  in  the  water),  provided  that  the  pH  is 
approximately  constant  over  the  study  area 
and  that  the  dissolved  [M]  in  the  water  to  which 
the  organisms  are  exposed  is  in  sorptive  equi- 
librium with  the  sediment  component  used  for 
normalization.  A  logical  consequence  is  also 
that  the  ratios  {M}  :  {Fe}  or  {M}  :  {OM},  al- 
though much  better  predictors  than  total  sed- 
iment metal  concentrations  in  some  cases,  are 
not  "universal"  predictors  in  the  sense  that 
their  application  should  be  site-dependent  (i.e. 
restncted  to  a  narrow  pH  range). 


Conclusions 

The  contamination  of  surficial  sediments  by 
metals  of  anthropogenic  origin  is  of  potential 
ecological  importance.  Past  attempts  to  relate 
total  metal  concentrations  in  surficial  sedi- 
ments to  metal  levels  in  indigenous  benthic 
organisms  have  generally  proven  disappoint- 
ing, as  noted  by  other  researchers  (Luoma  and 
Bryan  1982;  Bryan  1985;  NRCC  1988);  the 
present  study  offers  another  example  of  the 
failure  of  this  approach  (see  Fig.  4).  Even  when 
apparently  successful  within  a  limited  geo- 
graphical area,  such  empirical  models  cannot 
be  used  to  extrapolate  beyond  the  original  data 
used  to  develop  the  underlying  statistical  re- 
lationships. 

We  submit  that  a  more  promising  approach 


16 


Tessier  et  al. 


to  predicting  trace  metal  concentrations  in 
benthos  (and  eventually  to  predicting  metal 
effects  on  this  community)  would  be  to  de- 
velop deterministic  models  based  on  sound 
chemical  and  biological  principles.  Provided 
such  models  explicitly  incorporate  the  influ- 
ence of  major  environmental  variables  on  met- 
al uptake,  they  should  be  generally  applicable, 
i.e.  they  should  be  able  to  predict  the  metal 
concentration  in  benthos  from  lakes  other  than 
those  used  in  calibration.  The  development  of 
this  type  of  model  should  be  a  prerequisite  for 
use  of  animals  as  bioindicators  of  trace  metal 
pollution.  Such  an  approach  has  been  hindered 
in  the  past  by  our  incomplete  knowledge  of  the 
geochemical  and  biological  processes  that  con- 
trol metal  accumulation  in  aquatic  organisms. 

The  present  study  represents  an  attempt  to 
develop  such  a  deterministic  model  of  general 
applicability  to  predict  the  concentrations  of 
Cd  in  a  freshwater  bivalve.  In  formulating  the 
model,  we  used  concepts  derived  from  the  free- 
ion  activity  model  of  metal-organism  inter- 
actions and  from  surface  complexation  theory 
to  relate  Cd  concentrations  in  the  bivalve  to 
those  in  the  water  or  in  the  surficial  oxic  sed- 
iments. From  a  practical  point  of  view,  this 
approach  has  yielded  predictive  equations  for 
Cd  bioaccumulation  in  A.  grandis  that  require 
as  input  geochemical  variables  readily  mea- 
sured in  bottom  waters  and  sediments  (Table 
4).  The  predictive  power  of  these  equations  is 
acceptable  within  the  350,000-km^  study  area, 
but  they  should  obviously  be  tested  against  Cd 
accumulation  in  A.  grandis  in  lakes  in  other 
geographical  areas  before  their  generality  can 
be  proclaimed.  In  principle,  this  type  of  model 
should  apply  to  other  trace  metals  and  to  other 
sedentary  aquatic  organisms. 

Refinement  of  the  proposed  model,  or  for- 
mulation of  better  ones,  will  require  better  un- 
derstanding of  the  geochemical  and  biological 
processes  involved  in  bioaccumulation  of  trace 
metals  by  benthic  organisms.  For  example,  the 
nature  of  the  "sorption"  reactions  that  are  re- 
sponsible for  the  presence  of  trace  metals  in 
oxic  sediments  is  still  poorly  understood;  the 
identity  of  the  solid  phases  responsible  for 
binding  trace  metals  is  also  largely  unknown. 
Similarly,  improved  knowledge  of  the  feeding 
strategies  and  physiology  of  benthic  animals 
would  aid  in  the  development  of  future  mod- 
els. 
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